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ABSTRACT 

The  wide-band,  echolocating  sonar  signals  of  bats  are  investigated  with  reference 
to  statistical  estimation  theory  to  ascertain  the  implications  of  this  class  of  signals. 
The  work  is  divided  into  four  major  sections:  (1)  The  concept  of  an  ideal  receiver 
for  the  simultaneous  measurement  of  range  and  velocity  of  a  target  is  extended  to 
cover  wide-band  signals;  (2)  angular  localization  of  a  target  by  means  of  a  fixed, 
gain-dispersive  antenna  in  the  presence  of  additive  Gaussian  noise  is  studied; 
(3)  an  ambiguity  diagram  computer  is  designed  and  built  which  calculates  the  sig¬ 
nal  function,  defined  in  Chapter  2,  of  an  arbitrary  waveform;  and  (4)  ambiguity 
diagrams  are  calculated  for  several  types  of  bat  sonar  signals  emitted  by  the  bat 
Myotis  lucifugus  and  one  cruising  signal  emitted  by  the  bat  Lasiurus  borealis. 
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Little  brown  bat  Myotis  lucifugus  catching  tossed  mealworm. 


ECHOLOCATION  WITH  WIDE-BAND  WAVEFORMS: 
BAT  SONAR  SIGNALS 


1.  INTRODUCTION 


1.1  SCOPE  OF  THE  STUDY 

In  recent  years,  there  has  been  much  speculation  about  the  effectiveness  of  the  sonar  system 
of  bats.  Most  of  the  work  on  bats  has  pertained  to  measuring  the  capabilities  of  the  bat  in  differ¬ 
ent  problem  situations  such  as  forcing  the  bat  to  avoid  wires  in  a  noise  field,  measuring  the  bat's 
abilities  to  catch  flying  or  thrown  targets,  and  measuring  evoked  responses  in  the  bat's  auditory 
nervous  system.  Some  attempts  have  been  made  to  measure  and  to  understand  the  significance 
of  the  various  echolocating  signals  used  by  various  species  of  bats. 

This  research  investigates  the  way  in  which  the  characteristics  of  the  narrow-band  ideal 
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receiver  theory,  worked  out  by  Woodward  and  others,  apply  to  the  case  of  octave-wide-band 
bat  signals.  A  wide-band  signal  function  that  characterizes  the  operation  of  an  ideal  receiver 
is  defined  and  the  properties  of  this  signal  function  are  studied.  The  problem  of  elevation  esti¬ 
mation  with  a  wide-band  signal  is  studied  from  an  ideal  receiver  point  of  view  and  two  simple 
examples  are  worked  out  to  show  the  capabilities  of  this  type  of  system. 

An  analog  computer  is  developed  to  calculate  the  ambiguity  diagrams  of  several,  typical, 
bat  echolocating  signals.  Signals  of  bats  are  obtained  and  the  measurement  capabilities  of  these 
signals  are  deduced. 

1.2  BACKGROUND  OF  THE  PROBLEM 

In  the  18th  century,  Lazaro  Spallanzani  demonstrated  that  bats  need  both  their  ears  and 
mouths  unimpaired  to  fly  without  bumping  into  things,  while  removal  of  a  bat's  eye  does  not  affect 
his  agility  in  avoiding  obstacles.  At  that  time  this  curious  behavior  "contradicted  common  sense 
and  reasonableness  based  on  normal  experience."  In  1809,  Montagu  wrote,  "...  it  might  fairly 
be  asked,  'Since  bats  see  with  their  ears,  do  they  hear  with  their  eyes?'  "  (see  Ref.  3,  pp.  62-64). 

In  1938,  Donald  R.  Griffin,  who  is  now  Professor  of  Zoology  at  Harvard  University  but  was 
then  an  undergraduate  there,  entered  the  physics  laboratory  of  George  W.  Pierce  carrying  a  cage 
of  bats.  Pierce  had  just  perfected  a  microphone  sensitive  to  ultrasound.  When  they  turned  it 
toward  the  cage  of  bats,  Pierce's  loudspeaker  emitted  clicks,  demonstrating  that  bats  produce 
ultrasound  (see  Ref.  3,  p.67). 

Later  experiments,  conducted  by  Griffin  and  Robert  Galambos,  demonstrated  that  the  bat 
was  indeed  using  a  sonar  similar  in  principle  to  the  radars  and  sonars  that  were  being  developed 
by  man  (see  Ref.  3,  p.  69). 

In  1958,  Griffin  reported  some  of  his  work  in  Scientific  American,^  claiming  unusually  high 
"efficiencies"  for  bats  compared  to  man-made  radars.  Although  there  were  some  basic  errors 
made  in  interpreting  the  data,  the  bat  still  seemed  to  do  very  well. 
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M.I.T.  Lincoln  Laboratory  became  interested  in  studying  the  bat's  ability  to  avoid  wires  in 
the  presence  of  loud,  Gaussian  jamming  noise.  Improved  experiments  by  Griffin,  McCue,  and 
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Grinnell  and  more  comprehensive  theory  demonstrate  that  the  bat's  ability  to  avoid  wires  is 

very  good  but  does  not  violate  the  limits  imposed  by  information  theory. 

There  are  many  different  species  of  bats  using  echolocating  sonar  systems.  It  is  presumed 

that  since  the  bats  have  had  literally  millions  of  years  (see  Ref.  3,  p.  6)  in  which  to  evolve  their 

sonar  that  they  would  have  developed  some  sort  of  an  optimum  system.  The  numerous  and  widely 

distributed  Vespertilionidae  use  a  type  of  pulsed  FM  system;  the  Rhinolophidae.  the  "horseshoe 

bats"  of  Europe,  have  developed  a  pulsed  CW  system.  While  the  Vespertilionidae  keep  their  ears 
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stationary  during  flight,  the  Rhinolophidae  move  their  ears  at  rates  up  to  30cps. 

Each  individual  species  has  its  own  distinguishing  characteristics  and  individual  diet.  The 
vampire  bat  (family  Desmodontidae) ,  found  only  in  tropical  South  and  Central  America,  exists 
on  a  diet  consisting  entirely  of  blood  obtained  from  large  and  relatively  stationary  prey.  The 
Vespertilionidae,  however,  have  the  problem  of  detecting  and  pursuing  the  small  and  elusive  in¬ 
sects  on  which  they  feed.  For  the  purposes  of  this  study,  the  little  brown  bat  (Myotis  lucifugus) 
has  been  trained  to  catch  thrown  mealworms.  A  photograph  showing  the  Myotis  catching  a  meal¬ 
worm  appears  as  the  frontispiece.  The  Rhinolophidae  also  pursue  and  catch  small  insects  and 
moths.  A  photograph  taken  by  F.A.  Webster  and  D.  R.  Griffin  of  a  greater  horseshoe  bat 
(Rhinolophus  ferrum-equinum)  in  the  act  of  catching  a  moth  is  shown  in  Fig.  1.  There  are  bats 
that  feed  on  the  nectar  of  flowers  and  bats  that  feed  on  fruit.  One  of  the  most  unusual  diets  for 
a  bat  is  that  of  the  fishing  bat,  Noctilio  leporinus.  The  Noctilio  makes  his  living  by  catching 


Fig.  1.  Horseshoe  bat  (Rhinolophus  ferrum-equinum)  capturing  moth  thrown  by  D.  R.  Griffin. 
(Photograph  is  reproduced  with  permission  of  F.A.  Webster.) 
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small  fish  that  are  swimming  near  the  surface  of  the  ocean.  A  photograph  of  a  Noctilio  in  the 
act  of  catching  a  piece  of  fish  is  shown  in  Fig.  2. 

It  is  supposed  that  the  sonar  system  of  each  species  of  bat  has  developed  to  some  sort  of  an 
optimum  for  the  particular  task  which  that  species  performs.  Since  the  echolocation  problems 
facing  a  vampire  bat  are  much  simpler  than  the  problems  facing  an  insectivorous  bat,  the  sonar 

system  of  the  latter  is  expected  to  be  more  complex  than  that  of  the  former. 
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P.  M.  Woodward  and  others  have  shown  that  the  best  that  an  echolocation  system  can  do 
in  the  way  of  detection  of  a  signal  in  additive  Gaussian  noise  is  to  pass  the  received  signal  through 
a  matched  filter.  The  time-delay  and  velocity  measurement  capabilities  of  such  an  "ideal  re¬ 
ceiver"  are  represented  by  the  signal  function  which  is  derived  from  the  transmitted  signal. 

The  signal  function  of  the  bat's  sonar  pulse  contains  the  relevant  information  about  the  range 
and  range-rate  resolution  of  his  signal  when  it  is  used  in  an  additive  background  of  white  Gaussian 
noise.  When  it  is  clear  what  type  of  signals  are  used  by  the  different  species  of  bats  for  the 
different  operations  that  they  perform,  some  clues  as  to  the  way  in  which  we  can  improve  our 
modern  radar  systems  may  appear. 

The  bat  has  a  more  difficult  problem  than  that  of  determining  the  range  of  his  target.  He 
must  simultaneously  determine  the  range,  the  azimuth  and  the  elevation  as  well  as  determine 
the  type  of  target.  Very  little  is  known  about  the  process  involved  in  determining  the  elevation 
of  the  target.  Although  it  is  hypothesized  that  the  bat  determines  range  from  the  target -echo 
delay  and  that  he  determines  azimuth  by  the  same  methods  that  are  used  by  humans  to  determine 
azimuth,  it  is  not  clear  how  the  bat  or  the  human  determines  elevation  of  a  sound  source.  Some 
believe  that  very  small  movements  of  the  head  are  used  to  obtain  two  base  lines  for  interaural 


Fig.  2.  Fishing  bat  (Noctilio  leporinus)  catching  piece  of  fish  from  indoor  pool  of  water. 
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time  delay  measurements.  From  these  two  interaural  delay  measurements,  the  elevation  of  the 
source  is  supposed  to  be  predicted.  This  does  not  seem  reasonable  for  humans  or  for  bats.  A 
human  can  determine  the  elevation  of  a  transient  noise  with  his  head  stationary  and  a  vespertil- 
ionid  bat  does  not  wiggle  his  head  during  the  pursuit  of  a  target. 
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Azimuthal  angular  location  by  humans  has  been  reported  by  Jongkees  and  von  der  Veer, 
Pierce,^  and  Mills Mills  has  shown  that  humans  can  resolve,  under  rather  ideal  conditions, 
interaural  time  differences  on  the  order  of  1  psec.  Very  little  work  has  been  reported  on  the 
abilities  of  humans  to  locate  the  elevation  of  a  sound  source.  In  the  work  reported  by  Jongkees, 
et  al.,  angular  location  by  unilaterally  deaf  subjects  was  investigated.  Although  he  does  not  ac¬ 
tually  mention  the  possibility  that  gain-dispersive  transducers  are  the  mechanism  of  the  angular 
location  scheme,  he  reported  that  distortion  of  the  shape  of  the  outer  ear  affected  the  ability  of 
the  subject  to  locate  objects.  It  is  reasonable  to  assume  that  the  gain  pattern  of  a  person's  ear 
changes  greatly  with  frequency.  It  is  therefore  plausible  that  a  person  could  use  this  variation 
in  gain  pattern  to  determine  the  angular  elevation  of  a  sound  source.  It  is  also  plausible  that  a 
bat  would  use  this  type  of  system  to  determine  the  elevation  of  his  targets.  An  ideal  model  of 
this  type  of  system  is  investigated  in  this  report. 

1.3  OUTLINE  OF  THE  REPORT 

The  report  is  divided  into  three  parts:  theoretical  derivations  (Chapters  2  and  3),  descrip¬ 
tion  of  the  electronic  equipment  and  techniques  that  were  used  (Chapter  4),  and  the  results  of  the 
study  (Chapter  5) . 

Chapter  2  contains  the  theoretical  basis  of  the  ambiguity  diagram  concept.  The  Woodward 
concept  of  the  ideal  receiver  for  range  and  range-rate  detection  is  extended  to  treat  wide-band 
signals.  The  signal  function  for  a  wide-band  bat  signal  is  developed  and  the  resolution  limitations 
of  the  system  are  discussed.  A  different  way  of  looking  at  the  ambiguity  diagram  of  a  pulsed 
FM  signal  is  derived  and  then  related  to  the  bat  sonar  problem. 

A  derivation  of  the  ideal-receiver  concept  as  related  to  angular  location  by  means  of  a  gain- 
dispersive  antenna  is  made  in  Chapter  3.  A  gain-dispersive  antenna  is  an  antenna  with  an  angular 
gain  pattern  that  changes  as  a  function  of  frequency.  A  signal  function  for  the  angle-locating 
system  is  defined  which  relates  the  uncertainties  in  time-delay  (range)  measurements  to  the  un¬ 
certainties  in  angular  measurements.  The  special  case  of  "window -function  antennas"  is  treated. 
A  window-function  antenna  is  an  antenna  whose  impulse  response  at  a  particular  angle  is  related 
to  the  "window  function"  of  the  antenna  by  a  proportionality  factor  and  a  compression  or  expan¬ 
sion  of  the  time  axis.  The  relationship  between  the  window  function  of  the  antenna  and  the  signal 
function  is  derived  and  two  special  examples  are  treated.  It  is  presumed  that  bats  might  use 
this  type  of  system  for  angular  determination,  but  no  definite  proof  has  been  established. 

An  analog  ambiguity  diagram  computer,  which  was  designed  to  calculate  the  wide-band  sig¬ 
nal  function  developed  in  Chapter  2,  is  described  in  Chapter  4.  This  computer  calculates  the 
ambiguity  diagram  of  an  arbitrary  signal  that  is  recorded  on  magnetic  tape.  The  signal  is  re¬ 
recorded  from  the  magnetic  tape  onto  a  magnetic  drum  that  has  two  signal  tracks,  one  for  the 
direct  signal  and  one  for  the  Doppler-shifted  signal.  The  two  signals  are  shifted  with  respect  to 
each  other  by  means  of  a  movable  playback  head  on  the  magnetic  drum.  The  signals  are  multi¬ 
plied  and  integrated;  then  the  output  is  displayed. 
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The  results  obtained  by  analyzing  several  different  bat  signals  with  the  ambiguity  diagram 
computer  are  presented  in  Chapter  5.  Photographs  of  the  actual  bat  signals  as  well  as  the  am¬ 
biguity  diagrams  are  presented.  The  range  and  the  range-rate  resolution  of  the  different  bat 
signals  are  presented  as  a  result  of  the  ambiguity  diagrams.  The  results  and  conclusions  of 
this  research  are  summarized  in  Chapter  6. 
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2.  WIDE-BAND  ECHOLOCATION  THEORY 


2.1  INTRODUCTION 

If  a  suitable  (and  simple)  model  is  made  of  the  sonar  system  of  the  bat,  one  can  ask  what  is 
the  best  that  an  ’’ideal  receiver”  can  do  with  the  bat's  signals.  Some  of  the  interesting  charac¬ 
teristics  of  ideal  receivers  that  have  been  developed  in  the  radar  literature  were  developed  by 
assuming  that  the  signals  are  the  narrow-band  signals  usually  encountered  in  radar  systems. 
Clearly  there  is  some  doubt  as  to  whether  or  not  these  characteristics  apply  for  octave-wide¬ 
band  bat  signals.  Most  of  the  basic  derivations  in  the  literature  are  general  in  nature  and  apply 
to  all  types  of  signals,  or  require  only  slight  interpretation  for  application  to  wide-band  signals. 
In  particular,  for  narrow-band  signals,  a  Doppler  shift  is  considered  only  as  a  change  in  center 
frequency.  For  wide-band  signals  a  Doppler  shift  must  be  considered  as  an  expansion  or  com¬ 
pression  of  the  time  axis  of  the  transmitted  signal. 

A  complete  derivation  is  given  in  the  following  section  to  introduce  a  wide-band  signal  func¬ 
tion  and  investigate  its  characteristics. 

2.2  MODEL  AND  DERIVATION 

The  model  chosen  for  this  analysis  consists  of  a  moving  target  at  some  distance  from  the 
receiver,  as  shown  in  Fig.  3.  The  signal  at  the  receiver  input  is 


y(t)  \Ta  s  [ of ( t  —  r ) ]  +  n(t) 


(2-1) 


where 


y(t)  =  received  signal 


s(t)  =  transmitted  or  expected  signal 


n(t)  =  additive  Gaussian  noise 
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time  delay  associated  with  the  range  of  the  target 


a 


Doppler  factor 


v  target  velocity  or  range  rate 
c  =  speed  of  propagation 


] 3~45~7480  1 


ANTENNA 


p(a,r|y(t))  - 


IDEAL 

RECEIVER 


y  (t ) 


TARGET 


r 


2 


Fig.  3.  Model  of  echolocation  system  for  determination  of  time  delay  and  Doppler  shift. 
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The  term  multiplying  the  signal  in  Eq.  (2-1)  makes  the  signal  term  have  an  energy  that  is 
independent  of  the  velocity  of  the  target.  This  is,  in  effect,  an  assumption  that  the  target  cross 
section  is  independent  of  the  velocity  of  the  target. 

There  are  many  factors  omitted  from  this  model:  amplitude  variation  of  the  received  signal, 
target  cross  section,  antenna  pattern,  etc.  These  factors,  although  present  in  any  real  echo- 
location  system,  do  not  add  to  the  clarity  of  the  ideal  receiver  derivation  and  are  therefore 
omitted. 

2.2.1  Ideal  Receiver 

An  ideal  receiver  is  defined  as  a  system  that  destroys  all  unwanted  information  in  a  received 
signal  and  retains  all  wanted  information.  In  the  system  described,  the  wanted  information  is 
the  range  (ct/2)  and  the  Doppler  factor  (a)  of  the  target,  all  other  information  is  unwanted.  Since 
the  additive  Gaussian  noise  is  statistical  in  nature,  the  output  of  the  ideal  receiver  is  also  of  a 
statistical  nature.  All  the  information  contained  in  y(t)  concerning  r  and  a  is  expressed  as  the 
joint  conditional  probability  density  function  for  r  and  a  given  y(t): 

P(t ,  a  |y(t))  .  (2-2) 

No  additional  amount  of  processing  of  the  received  signal  can  extract  any  more  information  con¬ 
cerning  the  values  of  t  and  a.  Also,  in  obtaining  this  probability  function,  there  has  been  no 
information  about  t  and  a  that  has  been  destroyed.  Therefore,  a  receiver  that  computes 
p( t ,  cv  |y(t)),  or  something  monotonically  related  to  this,  is  the  ideal  receiver  for  this  problem. 

The  computation  of  p(r,a  | y (t ) )  is  performed  by  using  a  probability  theory  identity 

p(a,b|c)  p(c)  =  p(a,b,c)  =  p(c|a,b)  p{a,b)  (2-3) 

p(a,  b  |  c)  =  p(c  1  a,  b)  .  (2-4) 

By  substitution, 

p(r,a|y(t))=  pfyfy  p(y(t)|r,a)  ,  (2-5) 

where  p(T,aO  is  the  probability  density  function  for  t  and  a  before  the  reception  of  the  received 
signal.  It  is  called  the  a  priori  probability  density  function.  The  probability  density  function 
p(y(t))  is  a  strange  function.  To  be  mathematically  rigorous,  y(t)  would  be  confined  to  a  finite 
bandwidth,  the  sampling  theorem  would  be  applied,  and  the  samples  of  y(t)  would  be  used  to  form 
a  vector  in  waveform  space.  The  function  p(y(t))  is  therefore  the  probability  per  unit  volume  of 
waveform  space  that  a  particular  y(t)  is  received.  It  is  a  scalar  function  with  a  vector  argument. 
Likewise,  p(y(t)  |  r,  a) ,  the  likelihood  function,  is  the  probability  per  unit  volume  of  waveform 
space  that  a  particular  y(t)  is  received  given  that  the  actual  target  was  at  t  moving  with  a  Doppler 
factor  a. 

2.2.2  Likelihood  Function 

Since  the  noise  is  added  to  the  received  signal,  and  the  form  of  the  signal  is  known  exactly, 
the  likelihood  function  is  equal  to  the  probability  density  for  the  corresponding  noise  signal 

p(y(t)  I T,  a)  =  p(n(t))  .  (2-6) 
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The  notation  used  here  is  slightly  ambiguous  since  the  probability  function  notation  is  used  to 
denote  two  different  functions.  What  is  implied  by  this  notation  is  that  the  probability  density 
(probability  per  unit  volume  of  waveform  space)  for  a  particular  received  signal  y(t)  is  equal  to 
the  probability  density  for  the  particular  noise  signal  n(t)  that  was  received.  As  developed  by 
Woodward  (see  Ref.  27,  p.  38),  the  probability  density  function  for  a  particular  white  Gaussian 
noise  waveform  with  mean  energy  Nq  watts/cps  (single -sided  spectrum)  is  given  by 


p(n(t))  =  k  exp 


n 


<t)  dt 


(2-7) 


The  value  of  the  constant  k  depends  on  the  limits  of  the  integration.  From  Eq.  (2-1), 
n(t)  =  y(t)  -  s  [a(t  -  t)] 

Combining  this  with  Eqs.  (2-5),  (2-6),  and  (2-7), 


(2-8) 


P(r,  a  |y(t))  =  k  exp  [- ^  y2(t)  dt 

X  exp  ^  y(t)  s  [a(t  -  t)]  dt  exp  J  s2  [a(t  -  r)]  dtj  .  (2-9) 

All  the  factors  not  containing  t  and  a  can  be  lumped  together  in  one  common  constant  k.  The 
signal  is  normalized  by  setting 


u(t)  =  ——  s(t) 

\Te 


where  E  is  the  energy  of  the  signal,  given  by 

-I 


s2(t)  dt 


Combining  (2-9),  (2-10),  and  (2-11), 


P(t,  a  |y(t))  =  kp(r,  a)  exp  ^-j  exp  J  y(t)  u  [a(t  -  t)]  dtj 


(2-10) 


(2-11) 


(2-12) 


2.2.3  Receiver  Operation 

The  operation  of  the  receiver  on  the  received  signal  is  contained  in  the  exponent  of  the  last 
factor  of  Eq.  (2-12).  This  exponent  is  called 

q(T,  a)  =  2^**  f  y(t)  u  [a(t  -  t)]  dt  .  (2-13) 

o 

From  Eq.  (2-1),  y(t)  is  the  sum  of  two  terms,  the  signal  term  and  the  noise  term.  Assuming 
that  the  target  was  actually  at  tq  with  a  Doppler  factor  of  aQi 

q(T,  a)  =  h(T,  a)  +  g(r,  a)  (2-14) 

h(r,  a)  =  f  n(t)  u  [a(t  -  r)]  dt  (2-15) 

o  ^ 
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g <t,  a)  = 


(2-16) 


2EnP% 


N 


J  u  [aQ(t  -  tq)]  u  [a(t  -  r) 


dt 


Since,  in  Eq.(2-15),  n(t)  is  white  Gaussian  noise,  h(r,a)  is  a  Gaussian  variable  where  the 
mean  =  0  and 


2  2 

variance  =  oy  =  E  [h  (T,a)] 


(2-17) 


2  4Ea 

h "  V2" 


E  [n(t)  n(t  +  77)]  u  [a(t  —  t)]  u  [a(t  4-  77  —  r)]  dt  dr]  .  (2-18) 


The  noise  process  is  stationary  so  that  the  ensemble  average  is  equal  to  the  autocovariance 
function  <p(r] ).  In  Helstrom  (see  Ref.  21,  p.50),  the  autocovariance  function  for  white  Gaussian 
noise  is  given  as 

N 

<PIV)  =  E  [n(t)  n(t  +  7i)]  =  -f  >  (2-19) 

where  6(77)  is  the  Dirac  6-function.  When  the  77  integration  is  performed  on  Eq.(2-18)  this  yields 


2  _  2Ea 
o 

2  2E 
o 


ah  =  N 


ah  =  N 


^  u  [a(t  —  t)]  u  [a(t  —  t)]  dt 


(2-20) 


(2-21) 


The  ratio  2E/Nq  keeps  showing  up  in  these  equations.  This  is  the  signal-energy-to-noise-power- 
density  ratio  and  occurs  often  enough  to  be  given  the  special  name 

2E 


R  = 


N 


(2-22) 


dt 


Summarizing  the  previous  results: 

q (t,  a)  =  h(T,  a)  +  g (t,  a) 

h(T,  a)  =  Gaussian  random  variable  (mean  =  0,  cr^  =  R) 
g(T,  a)  =  R  J  u  [aQ(t  -  to)J  u  [a(t  -  r)]  Jaa~ 

2.2.4  Signal  Function 

The  integral  of  Eq.  (2-24)  can  be  simplified  by  setting 

t 1  =  a  (t  —  t  ) 
o  o 

r'  -  a  (t  —  t  ) 
o'  o 


[(2-14)] 


(2-23) 


(2-24) 


(2-25) 


then 
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g(  r of  * )  R  \  u(t’)  u  [a  1  ( t 1  -  t  1 )  ]  \fa~'  dt 


(2-26) 


This  integral  contains  the  interesting  characteristics  of  the  detection  processes  that  are  a  func¬ 
tion  of  the  transmitted  signal.  The  primes  are  removed  from  the  variables  in  Eq.  (2-26)  to 
simplify  the  notation  and  the  integral  is  called  the  signal  function 

X(r,  a)  y  u(t)  u  [ o' (t  —  r ) ]  \Ta  dt  .  (2-27) 

The  significance  of  x(t.o')  can  be  seen  by  substituting  back  into  the  probability  expression  (2-12) 

p(r,  a  |  y  (t ) )  =  kp(  t  ,  cx)  e^T,Q^  e^X<T,o?)  .  (2-28) 

If  the  a  priori  probability  density  function  p(rJa')  is  flat  and  for  a  given  R,  the  signal  function 
X(r,a)  is  proportional  to  the  logarithm  of  p(T,  a  | y(t)) ,  except  for  the  amount  contributed  by  the 
function  h(T,«).  Since  the  average  value  of  h(r,Q')  is  zero,  the  greater  x(T,a)  is,  the  more 
probable  that  particular  t  and  a  are.  x('J'.Q')  is  a  cross-correlation  function  between  two  signals 
with  unit  energy.  Therefore, 

X(0,  1)  1  >  x(t,  a)  .  (2-29) 

The  proportionality  factor  between  x(T,a)  and  In  [p(r,  a?  |y(t))  ]  is  R.  This  means  that  at  very 
high  signal-to-noise  ratios  (R  1)  very  good  r  and  a  estimation  is  possible. 

An  alternate  method  of  determining  the  signal  function  can  be  derived  by  substituting  the 
Fourier  transform 


u(t)  J  U(w)  eja,t  ^  (2-30) 

for  the  signal  in  Eq.  (2-27), 

X(t,  a)  U(a) ')  e^'1  ^  ^  U(w)  eja,a'<t_r)  ^  4a  dt  (2-31) 

=  OT  U(a>')  U(co)  ej(a,'+aw)t  e'jWaT  ^  ^  ^  dt  (2-32) 


but 

^  ej(«'+aw)t  dt  =  27r6(w,  +  aw)  (2-33) 

so  that 


X(r,  a)  =  J  U(-aw)  U(w)  e~iau>T  \Ta  ^ 

but  since  u(t)  is  real: 

U(-aw)  =  U*  («o>) 

Combining  (2-34)  and  (2-35),  the  alternate  form  for  x(t ,  a)  is  obtained: 


(2-34) 


(2-35) 
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X(t>  a)  =  J  U*(aoj)  U(a>)  e‘^awT  sfa 


(2-36) 


2.3  RANGE  AND  RANGE-RATE  RESOLUTION  LIMITATIONS 

OF  THE  A  POSTERIORI  MODEL 

There  are  several  conditions  which  limit  the  type  of  function  that  x(T,a)  can  be.  Since  this 
function  in  fact  characterizes  the  range  and  range  rate  resolution  limitations  of  an  echolocation 
system,  the  more  that  can  be  said  about  the  function  the  better  the  designer’s  position  is  when 
selecting  a  transmitted  signal  waveform. 

2.3.1  Probability  Limitation 

Since  p(a,r|y(t))  is  a  probability  density  function,  the  volume  of  this  function  must  integrate 
to  one.  That  is. 


(2-37) 


Referring  to  Eq.  (2-28)  it  is  seen  that  this  condition  tells  us  nothing  about  the  a  posteriori  model 
but  is  merely  a  way  of  evaluating  the  constant  k. 

2.3.2  Signal  Function  Area 

As  was  noted  in  Eq.  (2-29),  the  signed  function  for  r  =  0  and  a  =  1  is  equal  to  one.  Since  the 
height  of  the  signal  function  is  related  to  the  likelihood  that  the  target  is  at  some  range  traveling 
at  some  velocity,  if  the  signal  function  is  equal  to  one  at  some  point  other  than  (0,  1),  this  point 
in  range  and  range  rate  space  is  as  likely  as  the  point  (0,1).  In  other  words,  over  the  region 
that  the  signal  function  has  a  value  close  to  one,  the  values  of  range  and  range  rate  are  indis¬ 
tinguishable  or  ambiguous.  Therefore,  the  most  desirable  echolocation  signal  would  be  one  for 
which  the  signal  function  takes  on  a  value  of  one  at  (0,1),  and  zero  or  a  negative  value  elsewhere. 
Such  a  signal  is  not  possible  because  of  the  restrictions  on  the  form  of  x(t,  a). 

An  estimate  of  the  spread  of  the  signal  function  in  frequency  and  time  can  be  made  by 
calculating 


oO 


(2-38) 


In  order  to  calculate  this  function,  several  constants  are  defined.  If  we  assume  that  our  echo- 
location  signal  does  not  contain  any  DC  component. 


(2-39) 


and  has  unit  energy 


(2-40) 


a  center  frequency  o>o  can  be  defined  as  the  average  frequency  of  the  signal 
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(2-41) 


-o-r 


|  co  I  U*(oj)  U(w) 


do; 

2tt 


A  bandwidth  p  is  then  defined  as  twice  the  second  moment  of  the  signal  transform  around  the 
center  frequency 


(-§)2  =  r2=f  <|o,|  -co  )2  U*(«)  U(«)  .  (2-42) 

— OO 

The  absolute  value  signs  on  the  c j  are  needed  to  take  care  of  the  negative  half  of  the  signal 
transform. 

The  expression  (2-38)  can  be  rewritten  by  substituting,  from  Eq.  (2-36), 

oo 

U(au>)  U*(oj)  e^awT 


x[j~  U*(aaj ')  U(w')  e-jaw'T  ^  ^  dr 

Performing  the  t  integration, 

C  w  adr  =  2tt6(cj  —  cj!) 

J  — oo 

The  cj'  integration  can  also  be  performed,  reducing  (2-43)  to 

OO 

rr  h  00  ^ a 

U(aw)  U*(w)  U*(aco)  U(«)  -§j-  . 

— oo 


(2-43) 


(2-44) 


(2-45) 


Let 


then 


V  =  CJC* 


dv  =  coda 


(2-46) 


(2-47) 


The  integration  limits  on  v  are  the  same  as  on  a,  if  w  >  0.  If  cj  <  0,  the  integration  limits 
in  Eq.  (2-45)  are  reversed.  Changing  the  limits  back,  Eq.  (2-47)  is  written  . 


dy  =  —coda  (cj  <  0) 

Equivalently, 

dy  =  |  cj  |  da 

Equation  (2-45)  then  becomes 


U*(«)  11(0,)  U*00  L»  I7  ,4  ^  • 


(2-48) 


(2-49) 


(2-50) 


If  U(cj)  is  nonzero  only  near  cj  =  ±ojo,  cjq/|cj  |  can  be  expanded  in  a  Taylor  series  around 

both  +  cj  and  —  cj 
o  o 
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(cj  >  0) 


o 

M 


1 - (cj  —  CJ  )  +  — »  (oj  —  cj  )  — .  .  . 

cj  o  2  o 

O  CJ 

o 


1  +  -  (cj  +  CJ  )  -f  - j  (cj  —  CJ  )  +  .  .  .  (cj  <  0) 

CJ  O  £  o 

O  CJ 

o 


Reducing  this  to  a  simpler  form. 


cj  -  .  ~ 

T~T  =  1  — -  ( I  cj  I  —  cj  )  +  — 5  ( I  cj  I  —  cj  )  — ...  , 

CJ  '  1  1  n'  7  '  1  1  #V  1 


o'  2 
^o 


substituting  this  into  Eq.(2-50)  and  applying  Eq.  (2-40)  on  the  v  integration. 


(2-51) 


(2-52) 


f  U*(w)  U(w) 

J  -oo 

1  -7T  "“o'  +-S  (|a>l  "  wo)2  "• 

0 

3 

0 

dcj 

27T 


(2-53) 


By  substituting  Eqs.(2-40),  (2-41),  and  (2-42)  into  (2-53),  the  result  is 

1 2 


lx(T,ar)| 


9  ^  da 

2  O 

27T 


dr  =  1  + 


1C. 

. 


(2-54) 


This  equation  has  been  derived  by  integrating  a  from  to  .  For  a  less  than  zero,  the 
target  is  moving  toward  the  transmitter  at  a  velocity  greater  than  the  velocity  of  propagation  in 
the  medium.  For  sonar,  the  target  would  be  moving  toward  the  transmitter  at  a  velocity  greater 
than  the  velocity  of  sound.  For  radar,  the  target  would  be  moving  toward  the  transmitter  at  a 
velocity  greater  than  the  speed  of  light! 

Clearly  the  a  integration  should  be  performed  only  over,  at  most,  the  region  where  a  lies 
between  0  and  °°  .  This  restriction  is  easy  to  include  by  referring  to  Eqs.  (2-50),  (2-40),  and 
(2-46).  Since  u(t)  is  a  real  function, 

U*(i>)  U(J')  =  U(-v)  U*(-i>)  .  (2-55) 


The  argument  of  the  v  integration  in  Eq.  (2-50)  is  therefore  symmetrical  around  v  -  0.  The 
integral  on  a  between  0  and  <»  is  therefore  equal  to  one-half  the  integral  on  a  between  — «  and 
+  <»  .  The  result  as  stated  in  Eq.  (2-54)  is  then  changed  to 


Ix(t,q’)|2 


cj  da 
o 

2tt 


(2-56) 


This  is  an  important  result  and  requires  some  comment.  Equation  (2-56)  implies  roughly 
that  the  volume  under  the  signal  function  is  a  constant.  Roughly  speaking,  different  transmitted 
signals  just  change  the  position  of  the  fixed  amount  of  volume  contained  by  the  signal  function. 
This  is  almost  analogous  to  moving  around  a  fixed  amount  of  sand  in  a  box.  The  analogy  is  not 
perfect  because  Eq.  (2-56)  is  actually  a  statement  about  the  volume  under  the  square  of  the  signal 
function,  rather  than  the  actual  volume  under  the  signal  function.  Negative  values  of  the  signal 
function  are  therefore  treated  the  same  as  the  positive  values.  In  the  actual  measurement,  the 
negative  values  of  the  signal  function  signify  a  lower  likelihood  than  would  a  value  of  zero  and 
should  not  be  treated  the  same  as  positive  values  of  the  signal  function.  With  narrow-band  sig¬ 
nals,  the  peaks  and  the  valleys  of  the  signal  function  are  very  close  together  and  because  of  other 
limitations  of  the  transmitter  and  receiver  are  often  indistinguishable. 
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One  method  that  is  very  useful  in  getting  rid  of  ambiguous  area  in  the  a,  r  plane  is  to  use 
a  system  for  which  the  ambiguous  area  is  located  outside  the  limits  imposed  by  the  a  priori 
probability  conditions.  For  example,  the  operation  of  a  system  that  has  a  priori  constraints  on 
the  velocity  of  its  targets  can  be  improved  by  locating  an  ambiguous  area  outside  these  maximum 
velocity  limits. 

2.3.3  Ambiguity  Diagram  of  Pulsed  FM  Signals 

The  general  form  of  the  ambiguity  diagram  can  be  predicted  from  the  general  form  of  some 
transmitted  signals.  A  long-duration  constant-frequency  pulse  has  an  area  of  ambiguity  that  lies 
along  the  time  delay  axis  of  the  ambiguity  diagram,  whereas  a  short-duration  constant-frequency 
pulse  has  an  area  of  ambiguity  that  lies  along  the  Doppler  axis  of  the  ambiguity  diagram.  This 
meaas  that  with  the  former  pulse,  good  Doppler  resolution  is  achieved  at  the  expense  of  time 
resolution  and  with  the  latter  pulse,  good  time  resolution  is  achieved  at  the  expense  of  Doppler 
resolution.  With  a  pulse  that  is  long  in  duration,  with  the  frequency  varying  linearly  as  a  function 
of  time,  the  area  of  ambiguity  lies  along  a  line  that  makes  an  angle  with  the  time-delay  axis  of 
the  ambiguity  diagram.  As  stated  previously,  if  a  priori  constraints  on  the  echolocation  situa¬ 
tion  eliminate  the  possibility  that  the  target  is  moving  at  a  greater  velocity  than  some  fixed 
amount,  the  area  of  ambiguity  can  sometimes  be  reduced  by  substantial  factors.  Even  after  the 
a  priori  constraints  have  been  applied,  some  improvement  in  the  apparent  time -delay  measure¬ 
ment  accuracy  of  the  system  can  be  obtained.  Suppose  that  an  echolocation  system  that  uses  a 
pulsed  FM  signal  has  the  ambiguity  diagram  shown  in  Fig.  4  with  the  a  priori  constraints  also 
shown  in  the  figure.  The  basic  uncertainty  At  of  the  time  delay  estimation  of  this  system  is 
given  by  the  distance  between  the  two  dotted  lines.  This  deterioration  in  the  time  resolution  of 
the  system  over  the  resolution  that  is  obtained  when  the  target  velocity  is  known  exactly  is  caused 
by  asking  the  wrong  question  of  the  system.  If  the  question  were  asked  "Where  will  the  target  be 
in  t  seconds?"  rather  than  "Where  is  the  target  now?",  an  improvement  in  the  apparent  time 
resolution  of  the  system  can  be  obtained. 


Fig.  4.  Ambiguity  diagram  of  pulsed  FM  signal  with  a  priori 
velocity  constraints. 
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Let  the  slope  of  the  line  that  goes  through  the  center  of  the  area  of  ambiguity  in  Fig.  4  be 
p.  The  equation  for  this  line  is  then  given  by 

<*  -  1  =  jS(r  -  tq)  .  (2-57) 

Using  the  approximate  form  for  the  relation  between  the  target  velocity  v  and  the  Doppler  factor 
«  =  ,  (2-58) 


Eq.  (2-57)  can  be  reduced  to  an  expression  for  the  time  delay  associated  with  the  range  of  the 
target 


T  = 


(2-59) 


The  equation  of  motion  for  the  range  x  of  the  target  is 
X  =  X  +  vt 

o 


(2-60) 


where  xq  is  the  position  of  the  target  at  time  t  =  0.  The  position  of  the  target  at  time  t  ~  0  as 
determined  by  the  echolocation  system  from  Eq.  (2-59)  is 


cr 

cr  _ o  _  _v 

Xo  "  2  =  2  p 

Combining  Eqs.(2-60)  and  (2-6l), 


x  = 


+  vt 


which,  at  time 


yields  for  the  position  of  the  target: 


x  = 


(2-61) 


(2-62) 


(2-63) 


(2-64) 


Fig.  5.  Ambiguity  diagram  of  pulsed  FM  signal  modified  to  show 
ambiguity  diagram  for  position  of  target  at  later  time  (t  =  1/p). 
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independent  of  the  velocity  of  the  target.  Therefore,  the  position  of  the  target  can  be  determined 
as  accurately  as  if  the  velocity  of  the  target  were  known  exactly.  The  time  for  which  the  position 
of  the  target  is  known  is  a  little  later  or  a  little  earlier  than  the  time  that  the  signal  bounces  off 
the  target,  depending,  respectively,  on  whether  (3  is  positive  or  negative. 

The  modified  ambiguity  diagram  for  the  range  and  the  velocity  of  the  target  at  the  earlier 
time  is  shown  in  Fig.  5.  As  an  example,  a  radar  is  tracking  a  missile.  The  waveform  is  pulsed 
FM,  sweeping  up  in  frequency: 


Center  frequency  1  Gcps 

Frequency  sweep  10  Mcps 

Pulse  duration  10  psec 


Under  these  conditions,  a  1-msec  delay  will  cause  the  axis  of  the  ambiguous  area  to  rotate 
enough  to  coincide  with  the  r  t  zero  line. 
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3.  DIRECTIONAL  LOCATION  WITH  GAIN-DISPERSIVE  ANTENNAS 


3.1  INTRODUCTION 

A  bat  that  is  catching  a  target  needs  three  coordinates  of  information  about  his  target: 
range,  azimuth,  and  elevation.  He  can  obtain  range  information  from  the  target  echo  delay.  He 
can  obtain  azimuth  information  from  the  interaural  delay  of  the  echo.  The  problem  remains  as 
to  how  the  bat  obtains  elevation  information. 

One  theory  is  that  the  bat  wiggles  his  head  back  and  forth  to  obtain  several  interaural  delay 
measurements  with  different  baselines.  Although  the  bats  of  one  family,  the  Rhinolophidae,  are 
known  to  wiggle  their  ears  at  a  high  rate,  they  use  very  narrow-band  signals  which  would  give 
them  poor  time  resolution.  High-speed  photography  of  vespertilionid  bats,  in  particular,  the 
little  brown  bats  (Myotis  lucifugus)  in  pursuit  of  thrown  targets,  shows  no  such  head  wiggling, 
while  it  clearly  shows  that  the  bats  have  excellent  elevation  aim. 

Another  possible  explanation  is  that  the  bat  is  using  the  gain  pattern  of  his  ears  to  determine 
the  elevation  of  his  target.  The  bat's  ears  have  a  different  Fourier  system  function  for  different 
elevation  angles.  In  other  words,  the  gain  pattern  of  the  bat's  ears  varies  as  a  function  of  fre¬ 
quency.  The  bat  must  determine  from  what  angle  his  echo  is  returning  by  determining  which 
system  function  it  is  more  probable  that  his  signal  was  passed  through. 

The  problem  reduces  to  one  of  determining  the  impulse  response  of  a  linear  system  from 
the  noisy  output  of  the  system,  given  the  form  of  the  input  signal  and  the  class  of  functions  to  # 
which  the  impulse  response  belongs.  By  assuming  additive  Gaussian  noise  in  the  reception  sys¬ 
tem  that  is  independent  of  the  position  of  the  target,  the  problem  can  be  treated  in  the  same 
framework  as  the  problem  in  Chapter  2.  The  best  that  an  ideal  filter  can  do  is  to  compute  the 
a  posteriori  probability  distribution  for  the  elevation  angle  and  the  time  delay  of  the  target,  given 
the  received  signal.  An  antenna  function  is  defined  which  is  the  complex  Fourier  system  function 
for  the  antenna  at  the  angle  9.  A  signal  function  is  then  defined  which  describes  the  operation  of 
the  time  delay  and  elevation  measurement  system.  This  signal  function  is  similar  to  the  signal 
function  described  in  Chapter  2,  except  that  this  signal  function  depends  on  the  form  of  the  an¬ 
tenna  as  well  as  on  the  form  of  the  transmitted  signal. 

The  special  case  of  a  window-function  antenna  is  examined.  A  window-function  antenna  is 
an  antenna  which  has  a  sensitive  area  that  lies  on  a  plane  and  has  a  local  sensitivity  that  can  be 
defined  as  a  function  of  the  position  on  the  plane  of  the  antenna.  Any  one  point  on  the  antenna  is 
assumed  to  be  equally  sensitive  to  a  signal  from  any  direction. 

Two  window-function  antennas  are  used  as  examples  to  demonstrate  the  concept  of  the  an¬ 
tenna  function,  the  signal  function  and  the  operation  of  the  echolocation  system. 

3.2  MODEL  AND  DERIVATION 

The  model  used  for  this  analysis  is  shown  in  Fig.  6.  The  echo  signal  comes  into  the  receiv¬ 
ing  antenna  at  an  angle  0.  After  the  signal  is  received  by  the  antenna,  white  Gaussian  noise  is 
added.  Therefore,  the  received  signal  y(t)  consists  of  two  parts:  the  transmitted  signal  s(t) 
passed  through  a  linear  filter  (the  antenna),  and  the  additive,  white  Gaussian  noise 


-oo 


(3-1) 
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Fig.  6.  Gain-dispersive  antenna  model. 


where 


s(t)  =  transmitted  signal 

n(t)  =  additive  Gaussian  noise 

t  =  time  delay  associated  with  the  range 
h(0,t)  =  impulse  response  of  antenna  at  angle  0 
y(t)  =  received  signal 

If  the  Fourier  transform  of  s ( t)  and  h(0,  t)  are  used. 


s(t) 

jj  S(cj) 

jo>t  da) 

(3-2) 

h(e,  t)  = 

f  H(0. 

co )  e^1 

da> 

27T  ' 

(3-3) 

Eq.  (3-1)  can  be  written  in  the  form 

y(t)  =  n(t)  +  S(u)  ej“<r?"T)  £  H(0,  w')  eja,,(t'I,)  ^  drj  .  (3-4) 


By  noting  that 


C  eJ(^  <*>  ^  d?7  =  27r6(cj  —  u') 

_oo 

and  performing  the  a>'  integration,  Eq.  (3-4)  reduces  to 
y(t)  =  n(t)  +y  S(w>  H(G.«)  eja,(t"T)  g 


(3-5) 


(3-6) 


3.2.1  Signal  Function 

The  a  posteriori  probability  density  function  for  0  and  r  is  given  by 

p(Q,  T  | y(t) )  =  p(y<t)l  e.  T)  ,  (3-7) 

where,  since  the  convolution  integral  in  Eq.  (3-1)  is  entirely  deterministic  when  the  parameters 
0  and  t  are  given,  the  probability  density  function  for  the  likelihood  function  is  given  by 


p(y(t)|o,  r)  =  p(n(t))  =  k  exp 


n2(t)  dt 


(3-8) 
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where  with  the  use  of  Eq.  (3-6)  the  exponent  of  Eq.  (3-8)  can  be  rewritten 


~W  J  dt  =  -fT 

n 


=  ~d~  [  y2(t)  dt  + 

n 


y  |y(t)  S(co)  H(e,«o)  eja;(t_T)  ^J2 
if-  Jy(t)  [J  s(«)  H(e,  u)  ejw(t"T)  ^ 


dt 


dt 


(3-9) 

(3-10) 


-jf-  S(co)  H(0,co)  S(co')  H(G,  co ')  ej(aj+w')(t'T)  ^  ^  dt  . 

The  first  term  in  Eq.  (3-10)  concerns  only  y(t)  and  can  therefore  be  lumped  with  p(y(t))  in 
Eq.  (3-7)  and  k  in  Eq.  (3-8).  The  third  term  is  related  to  the  energy  of  the  signal  from  the  an¬ 
tenna  that  is  expected  at  a  particular  angle.  Letting  this  term  be  called  £(0,  T)  and  after  per¬ 
forming  the  integration 

C  eJ ( c*-) )t  ^  _  2tt6(cj0  +  co 1 )  ,  (3-11) 

d  _oo 

the  integration  on  cj!  can  be  performed  leaving 

|(0,  t)  =-fr]  S{w)  H(e.w)  S(-co)  H(0,-co)  ^  •  (3-12) 

Since  both  s(t)  and  h(0,  t)  are  real 

S(—  o>)  =  S*(oj) 

H(0,-w)  =  H*(0,  w) 

Eq.  (3-12)  reduces  to 

S(e,  t)  =  s*(w)  s(w)  H*(e,w)  H(e,w)  ^  .  (3-15) 

Assuming  that  the  actual  target  is  at  the  angle  0q  and  at  the  range  corresponding  to  the  time 
delay  r  ,  the  second  term  of  Eq.  (3-10)  can  be  written  in  two  parts:  the  part  due  to  the  signal 
tf>(0 ,  t),  and  the  part  due  to  the  noise  r(0,  t) 

4>(e,  t)  =  jf  -  ^  [  j*H(eo,  co)  S(co)  ejc0<t  To)  ^ 

H(e,co')  S(co')  eJC0'^t-7^  dt 

•'(e.  T)  5  jf-  j  n(t)  [  J  H(0,  co)  S(co)  ejw(t'r)  |^]  dt 
The  t  integration  in  Eq.  ( 3  - 1 6)  is 

r  ej("'+co)t  dt  =  27r6(co'  +  co)  (3-18) 


(3-16) 

(3-17) 


(3-13) 

(3-14) 
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Performing  the  c o'  integration  and  using  Eqs.  (3-13)  and  (3-14),  Eq.  (3-16)  becomes 


4(0,  t) 


Z_ 

N 


o 


H*(e,  co)  H(eQ.  co) 


I  S(co)  |  2 


jw('-r 

e 


dc o 
Ztt 


(3-19) 


Since,  unlike  the  problem  of  Chapter  2,  the  expected  energy  of  the  received  signal  can  de¬ 
pend  on  the  parameter  that  is  being  measured  (0),  the  function  which  describes  the  operation  of 
the  ideal  receiver  on  the  signal,  the  signal  function  (q(0,  r)),  consists  of  the  sum  of  the  operations 
described  by  Eqs.  (3-15)  and  (3-19), 

q(0,  t)  =  4(0,  t)  4  4(0,  t)  (3-20) 

2  rl  ,  ,2  Jw(T‘To) 

q(e,  r)  =  |H*(6.w)  H(0q,  co)  |S(w)|  e 

-j  |H(0,w)|2  |S(w)|2]  ^  .  (3-21) 

Normalizing  the  signal 


S2(t)  dt  =  E 


u(t)  = 

1 

S(t) 

\Te 

UM  - 

1 

S(w) 

\Te 

R  = 

2E 

N 

o 

Eq.(3-21)  reduces  to 

q(0.  t)  H  ^  |h* 


,2  jw<T-To> 

(O,oj)  H(0o,co)  |  U(oj>  |  e 


1 

2 


H(0,  co)  |  2  |UM|2 


d co 
Ztt 


(3-22) 

(3-23) 

(3-24) 

(3-25) 


(3-26) 


This  q(0,  t)  function  is  similar  in  several  ways  to  the  x(^>  t)  function  of  Chapter  2.  The 
meaning  of  this  function  can  be  seen  by  substituting  Eq.  (3-26)  back  into  Eq.  (3-7)  by  using 
Eqs.  (3-8),  (3-10),  (3-15),  (3-16),  (3-17).  and  (3-20) 

p(0,  t |  y (t ) )  k  p(0,  t)  e eq^0,  .  (3-27) 

This  signal  function  q(0,  r)  is  proportional  to  the  logarithm  of  the  a  posteriori  probability 
density  function  for  0  and  t.  Since  the  logarithm  is  a  monotonically  increasing  function,  the 
probability  density  that  the  target  is  located  at  a  particular  ©  and  r  is  monotonically  related  to 
the  signal  function.  This  means  that  the  greater  the  value  of  the  signal  function,  the  greater  the 
corresponding  probability  density  function  will  be. 

The  form  of  the  signal  function  is  determined  by  the  antenna  function  and  the  transmitted 
signal.  Since  the  normalized  transmitted  signal  spectrum  (U(to))  appears  in  Eq.  (3-26)  only  as  the 
absolute  magnitude,  signals  with  identical  magnitude  of  signal  spectrum  will  produce  identical 
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signal  functions.  Therefore,  an  impulse,  white  noise,  and  swept  FM  signals  would  give  identical 
signal  functions.  This  is  not  to  say  that  if  the  receiver  is  set  for  use  with  an  impulse  as  the 
transmitted  signal,  white  noise  could  be  used.  The  receiver  has  to  be  set  up  for  the  exact  form 
of  the  transmitted  signal.  Under  these  conditions,  the  above  signals  will  give  identical  signal 
functions. 

3.2.2  Noise  Function 

The  noise  function  r(9,  r),  defined  by  Eq.  (3-17),  is  a  Gaussian  random  variable.  Since  n(t) 
has  zero  mean,  r(0,  r)  also  has  zero  mean.  The  variance  of  the  noise  function  is 

(Tn2  =  E  [VZ(Q,  T>]  =  E  (n(t)  n(t')  ] 

No 

x  H(e,cj)  H(e,w')  S(oj)  S(oj')  eJw(t"T)  e-’a,'(t'"T)  dt  dt ' 

but 

N 

E  [n(t)  n(t')  ]  =  -f  5(t  —  f)  , 
performing  the  t*  integration  and  then  the  integration  on  a;1 

^  ^  dt  =  Zn6(u  +  a ;') 

and  substituting  from  Eqs.(3-24)  and  (3-25),  Eq.(3-28)  reduces  to 

%  =  R  J  |H(e,a>)|2  |U(w)|2  ^  .  (3-31) 

3.3  WINDOW- FUNCTION  ANTENNAS 

The  class  of  antennas  that  can  be  described  by  window  functions  allows  some  simplification 
of  the  signal  function.  A  window-function  antenna  is  an  antenna  whose  sensitive  area  lies  on  a 
plane,  no  energy  storage  takes  place,  and  the  sensitivity  at  any  point  can  be  either  positive  or 
negative.  An  acoustical  example  of  this  type  of  antenna  is  the  piston-in-a-wall  transducer. 

3.3.1  Impulse  Response  of  a  One -Dimensional  Window- Function  Antenna 

The  point  sensitivity  of  the  antenna  is  g(x/c).  Each  point  on  the  antenna  is  acting  independ¬ 
ently  of  all  other  points  on  the  antenna  and  is  independent  of  frequency.  The  impulse  response 
of  this  antenna  at  an  angle  9  can  be  calculated  by  inspection  from  Fig.  7(a-b). 

The  relationship  between  a  point  of  time  on  the  impulse  response  and  the  coordinate  of  the 
window  function  is 

t  =  ^  sin  0  .  (3-32) 

Since  each  point  on  the  antenna  is  considered  separately,  the  area  under  the  impulse  response 
is  a  constant,  equal  except  for  a  sensitivity  constant  to  the  area  under  the  window  function. 
Choosing  the  scale  for  g(x/c)  in  such  a  way  that  the  sensitivity  constant  is  unity, 

h(e,t)  dt  =  g(f)  d(-f)  .  (3-33) 


(3-28) 


(3-29) 


(3-30) 


21 


Fig.  7(b).  Construction  showing  timing  of  impulse  coming 
from  angle  0  on  window-function  antenna. 

By  differentiating  Eq.  (3-32) 

dt  =  sin  0  d( — ) 
c 

and  substituting  into  Eq.  (3-33) 

h(0,  l)  =  sine  g*  slrTe  * 


(3-34) 


(3-35) 


which  is  the  relationship  between  the  impulse  response  of  a  window-function  antenna  at  a  par¬ 
ticular  angle  and  the  window  function.  Under  the  normalizing  condition  on  the  window  function 


y  g2(n)  dp.  =  1 

the  energy  of  the  impulse  response  is 


(3-36) 


s 


h  (0,  t)  dt  = 


.  2 

sin  0 


h  (0,  t)  dt  = 


sin  0 


f 


2. 

g  ( 


sin  0 


)  dt 


(3-37) 


(3-38) 


Notice  that  as  0  0,  the  energy  of  the  impulse  response  This  is  normal  and  expected, 

since  the  energy  in  a  unit  impulse  is  infinite. 
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3.3.2  Transform  Relationship 

The  Fourier  transforms  for  the  impulse  response  and  the  window  function  are 


H(0,w)  =  f  h(0,t)  e'-^*  dt 

J  -oo 

(3-39) 

p  00 

G(t?)  =  \  g(p.)  dfi 

J  -OO 

(3-40) 

Combining  Eqs.  (3-35)  and  (3-39) 

-y  Je  «( • 

(3-41) 

by  substituting  p  for  t/sin0  this  becomes 

H(0.U)  g(n)  e“J"  sinG^ 

(3-42) 

which,  with  Eq.  (3-40),  reduces  to 

H(0,  w)  =  G(cj  sinG)  .  (3-43) 


3.4  EXAMPLES  OF  THE  SIGNAL  AND  LIKELIHOOD  FUNCTIONS 

Two  examples  are  worked  out  to  she  «  graphically  the  possibilities  of  this  angular  measure¬ 
ment  scheme.  To  simplify  the  arithmetic,  the  transmitted  signal  spectrum  is  assumed  to  be  flat 

|S(w)|2=l  (3-44) 

This  assumption  gives  infinite  bandwidth  to  the  transmitted  signal,  obviously  not  the  case  in  any 
real  echolocation  system.  Because  of  the  nature  of  the  two  window-function  antennas  described 
below,  the  infinite  bandwidth  assumption  affects  the  signal  function  radically  only  when  Gq  is 
near  zero. 


3.4.1  Antenna  With  a  Square  Window  Function 

The  window  function  of  this  example  is  shown  in  Fig.  8.  The  equation  for  the  function  is 


8(f)  = 


1 

\l~a 

0 


(3-45) 


The  Fourier  transform  of  the  window  function  is 


G(tj)  =  \ fa 


<2£> 


From  Eq.  (3-43), 

H(0,  cj)  =  G(cj  sin  0) 


(3-46) 


[(3-43)] 
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Fig.  8.  Window  function  for  square-window- 
function  antenna. 


x/c 


Fig.  9.  Square-window-function  antenna:  impulse  response  for  target 

at  0  ,  t  and  for  target  at  0,t. 
o  o  ** 


Fig.  10.  Convolution  integral  for  impulse  responses 
of  square-window-function  antenna. 
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the  transform  of  the  impulse  response  is 


H(0,  c o)  =  'fa 


sin  ( co  sin  ©) 
( ^  co  sin  0) 


and  from  Eq.  (3-35), 


h(0,t)  sin©  g(sin©) 


the  impulse  response  is 


(3-47) 


[(3-35)) 


h(0,  t) 


1 

\fa  sin  0 
0 


t 


sin  0 


t 


sin  0 


I  > 


a_ 

2 


(3-48) 


Since  multiplication  in  the  frequency  domain  is  equivalent  to  convolution  in  the  time  domain, 
under  the  conditions  of  Eq.  (3-44),  Eq.  (3-21)  can  be  rewritten  as 

q(e,  T)  =  [h(©o,  p)  *  h(e,  P-(T  -  To»  -|-y  h2(©,t)  dt]  .  (3-49) 

From  Eq.  (3-38)  the  second  term  of  the  above  integral  is 


2 


s 


h  (0,  t)  dt  =  - 


1 


2  sin  © 


(3-50) 


The  antenna  impulse  responses,  as  shown  in  Fig.  9,  when  convolved  yield  the  function  shown  in 
Fig.  10.  A  three-dimensional  graph  of  the  signal  function  with  0q  =  30°  is  shown  in  Fig.  11. 

Sectioning  through  the  signal  function  with  the  plane  t  =  0,  the  equation  for  the  signal  function 
becomes 


q(©,0) 


2 

f  .  /  1 

1  \ 

1  1 

N 

0 

[mln(sin0o  ' 

sin  0 / 

2  sin  0  J 

(3-51) 


which  is  shown  graphically  in  Figs.  12  and  13  for  two  values  of  0q  along  with  the  corresponding 
likelihood  functions  for  various  signal-to-noise  ratios. 


3.4.2  Antenna  With  a  Cosine  Window  Function 

The  window  function  for  this  example  is  shown  in  Fig.  14.  The  equation  for  the  function  is 


g(f)  - 


fIcosl°l(X)  |i|  <  « 

n la  a  c  1  c 1  2 

o  ,  l^-l  >  4 


(3-52) 


and  from  Eq.  (3-35)  the  impulse  response  of  the  antenna  at  an  angle  0  is 


h(0,  t) 


J1  — 

v  a  sin 


7:  COS  10  7T  - : - - 

0  a  sin  0 


hi  < a  sine 
N  > 


2 

a  sin  0 


(3-53) 
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Fig.  11.  Three-dimensional  graph  of  signal  function  of  square-window-function  antenna. 


NORMALIZED  LIKELIHOOD  FUNCTIONS  SIGNAL  FUNCTION 


Fig.  12.  Signal  and  likelihood  functions  for  square- 
window-function  antenna  (t  =  0,  0.  =  30°). 


Fig.  13.  Signal  and  likelihood  functions  for  square 
window-function  antenna  (t  =  0,  0  =  60°). 
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Fig.  14.  Window  function  for  5-cycle-cosine 
window-function  antenna. 


Let 

/  1  1  \ 

a  max  gin0  •  sin  0  / 
'  o ' 

,  _  .  /  1  1  \ 

min  \sin  0  *  sin  0  / 

'  o ' 


(3-54) 

(3-55) 


Sectioning  through  the  signal  function  with  the  plane  r  =  0  the  equation  for  the  signal  function 
becomes 


q(e.o)  =  jp 


2  ab 


5?r(a  -  b  ) 


t  sin  57r - ^ — : — 7T 

2,  a  2  sin  0 


(3-56) 


which  is  shown  graphically  in  Figs.  15  and  16  for  two  values  of  0q  along  with  the  corresponding 
likelihood  functions  for  various  signal-to-noise  ratios. 


3.5  FINITE  BANDWIDTH  CONSIDERATIONS 

The  two  previous  examples  (Secs.  3.4.1  and  3.4.2)  were  worked  out  under  the  assumption 
that  the  transmitted  signal  had  an  infinite  bandwidth.  This  can  obviously  never  be  the  case  in 
an  actual  echolocation  system.  Equation  (3-49)  was  derived  from  Eq.  (3-26)  under  the  infinite 
bandwidth  assumption.  For  a  finite  bandwidth,  the  first  term  of  Eq.  (3-49),  the  cross-correlation 
function  of  the  antenna  impulse  response,  should  be  cross  correlated  with  the  normalized,  signal 
autocorrelation  function.  The  second  term  can  be  written  as  it  appears  in  Eq.  (3-26).  The  effect 
of  the  cross  correlation  with  the  autocorrelation  function  can  be  seen  intuitively.  Suppose,  for 
instance,  that  the  signal  autocorrelation  function  is  very  narrow,  and  has  the  general  form  of  a 
Gauss  error  function.  Cross  correlating  the  narrow  Gauss  error  function  with  the  function  rep¬ 
resented  by  the  first  term  of  Eq.  (3-49)  has  the  effect  of  rounding  the  corners,  in  the  T-direction, 
and  generally  smoothing  the  shape  of  this  function.  This  rounding  of  the  sharp  corners  in  the 
T-direction  will  usually  round  sharp  corners  in  the  0 -direction  also.  The  flat  regions  are  not 
affected,  since  the  Gauss  error  function  integrates  to  some  constant  value  over  these  regions. 

A  signal  that  has  a  wider  autocorrelation  function,  and  therefore  a  narrower  bandwidth,  will 
cause  a  greater  rounding  of  the  corners  of  the  first  term  of  the  signal  function.  In  the  region 
where  this  autocorrelation  function  of  the  antenna  impulse  response  is  narrow,  near  0q  =  0,  the 
signal  autocorrelation  function  will  affect  the  signal  function  radically.  The  signal  function  takes 
on  the  rough  shape  of  the  signal  autocorrelation  function  rather  than  the  shape  of  the  autocorrela¬ 
tion  function  of  the  antenna  impulse  response. 

In  any  real  echolocation  system,  limitations  on  the  lower  frequency  cutoff  must  be  imposed 
as  well  as  the  upper  frequency  limitations.  In  particular,  it  is  unreasonable  to  assume  that  the 
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NORMALIZED  LIKELIHOOD  FUNCTIONS  SIGNAL  FUNCTION 


Fig.  15.  Signal  and  likelihood  functions 
for  5-cycle-cosine  window-function  an¬ 
tenna  (t  =  0,  0  =  30°). 
o 


Fig.  16.  Signal  and  likelihood  functions 
for  5-cycle-cosine  window-function  an¬ 
tenna  (t  =  0,  0q  =  60°). 
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echolocation  signal  has  a  DC  component.  If  the  signal  does  not  have  a  DC  component,  the  signal 
autocorrelation  function  does  not  have  a  DC  component.  The  effect  of  this  restriction  can  be  ob¬ 
served  by  referring  to  Fig.  11.  Slices  parallel  to  the  9  =  0  plane  are,  in  effect,  passed  through 
a  band-pass  filter  whose  impulse  response  is  the  autocorrelation  function  of  the  transmitted  sig¬ 
nal.  (Since  the  autocorrelation  function  is  nonzero  for  negative  time  as  well  as  positive  time, 
this  hypothetical  band-pass  filter  is  physically  unrealizable.)  The  DC  level  for  these  signals,  or 
slices,  can  be  determined  by  referring  to  Eq.  (3-26).  The  second  term  of  this  equation  is  the  one 
that  controls  the  DC  level  of  the  slices.  Changing  the  equation  from  the  frequency  domain  to  the 
time  domain  and  using  correlation  notation,  the  second  term  of  Eq.  (3-26)  can  be  written  as 

—  h(0,  t  —  77)  *  h(9,  |i  —  77)  *  u(p  —  /?)  *  u(— p)  ,  for  t  =  0  ,  (3-57) 


where  77,  p,  and  f$  are  dummy  variables.  The  DC  level  of  the  output  of  the  hypothetical  band¬ 
pass  filter  for  each  of  these  slices  is  therefore  equal  to  minus  one-half  the  output  of  this  filter 
at  r  =  0  with  the  autocorrelation  function  of  the  antenna  impulse  response  on  the  input. 

The  graph  in  Fig.  11  is  changed  by  the  finite  bandwidth  restriction  in  the  following  ways: 

The  level  of  the  graph  goes  below  the  DC  level  as  well  as  above  it. 

The  flat  region  in  the  vicinity  of  60°  to  90°  sags  slightly. 

The  sharp  edges  are  rounded  and  the  central  spike  is  smoothed. 

The  DC  level  of  the  graph  is  modified  to  equal  the  level  given  by  Eq.  (3-57). 


The  graph  in 

(i) 

The 

(2) 

The 

(3) 

The 

(4) 

The 

30 


4.  AMBIGUITY  DIAGRAM  COMPUTER 


4.1  INTRODUCTION 


An  ambiguity  diagram  computer  is  a  device  that  computes  the  signal  function  expressed  by 


Eq.  (2-27) 


[(2-27)] 


The  terms  "  ambiguity  diagram"  and  "signal  function"  are  used  almost  interchangeably  in  the 
radar  literature.  P.  M.  Woodward  first  used  the  term  signal  function  to  describe  the  narrow- 
band  approximation  to  the  function  mentioned  above.  Since  this  function  describes  the  area  of 
ambiguity  of  a  received  signal  waveform,  it  has  also  been  called  an  ambiguity  function.  To 
avoid  ambiguities,  the  above  function  will  be  called  the  signal  function  while  the  graph  of  this 
function  will  be  called  the  ambiguity  diagram. 

The  basic  requirements  of  an  ambiguity  diagram  computer  are: 

(a)  A  device  to  simulate  and/or  repeat  the  signal  u(t), 

(b)  A  method  of  producing  a  variable  time  delay, 

(c)  A  method  of  producing  a  Doppler  shift, 

(d)  A  multiplier, 

(e)  A  pulsed  integrator, 

(f)  A  method  of  displaying  the  output. 

There  are  many  ways  to  obtain  each  of  the  above  operations.  One  rather  obvious  solution 
would  be  to  compute  the  signal  function  on  a  digital  or  an  analog  computer.  Several  years  ago, 
a  student  in  the  M.I.T.  acoustics  laboratory,  Kenneth  Goff,  built  an  analog,  magnetic  drum,  time 
delay  unit.^  This  unit  records  one  or  two  signals  on  two  tracks  of  a  magnetic  drum.  One  of  the 
record  heads  is  movable  so  that  the  signal  being  recorded  on  one  track  can  be  advanced  or  re¬ 
tarded  with  respect  to  the  signal  being  recorded  on  the  other  track.  Because  this  time  delay 
unit  was  available,  it  was  decided  that  a  special-purpose,  analog  ambiguity  diagram  computer 
would  be  constructed  to  analyze  the  echolocation  signals  of  bats. 

The  Goff  delay  drum  was  originally  built  to  continuously  record,  reproduce,  and  erase  on 
the  two  signal  tracks.  With  only  minor  modifications,  the  unit  was  redesigned  to  accept  a  pulse 
signal.  The  pulse  signal  is  recorded  on  the  drum  and  remains  there.  Since  the  drum  is  rotating, 
the  pulse  signal  is  repeated  over  and  over  again.  The  signals  on  the  two  tracks  are  advanced  or 
retarded  with  respect  to  each  other  by  advancing  or  retarding  a  movable  playback  head.  When 
the  correlation  function  has  been  plotted,  the  signals  are  erased  and  the  drum  is  ready  to  accept 
the  next  pair  of  pulse  signals. 

The  signals  of  the  bats  are  recorded  on  a  high-speed,  instrumentation  tape  recorder.  By 
reproducing  the  tape  at  l/32  of  the  recorded  speed,  the  20-  to  120-kcps  frequency  of  the  bat  sig¬ 
nals  is  reduced  to  a  frequency  range  of  600  to  4000  cps,  which  is  audible  to  human  ears.  It  thus 
falls  well  within  the  frequency  band  for  which  the  Goff  time-delay  unit  is  designed  to  work. 

Doppler  shift  is  introduced  by  changing  the  speed  of  the  tape  reproducer  by  incremental 
amounts.  This  Doppler  shift  is  a  compression  or  expansion  of  the  time  axis  of  the  signal,  thus 
simulating  an  actual  Doppler  shift. 
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A  vacuum-tube  quarter-square  multiplier  is  used  for  the  multiplication  of  the  two  signals. 

29 

The  multiplier  was  originally  designed  and  built  as  a  Master's  thesis  project  by  the  author. 

The  pulsed  integrator  was  built  out  of  Philbrick  operational  amplifiers  with  the  gating,  sampling, 
and  resetting  operations  performed  by  electronically  driven  relay  circuitry. 

The  output  of  the  computer  is  displayed  on  a  Variplotter  X-Y  recorder.  When  a  three- 
dimensional  graph  output  is  desired,  tan  30°  times  the  X  deflection  of  the  recorder  is  subtracted 
from  the  Y  deflection  of  the  recorder.  The  zero  line  of  the  output  of  the  computer  is,  therefore, 
plotted  at  a  30°  angle  with  respect  to  the  horizontal.  A  tri-metric  drawing  is  then  traced  from 
a  composite  of  these  slanted  graphs  to  obtain  a  two-dimensional  representation  of  these  three- 
dimensional  graphs. 

4.2  RECORD  ELECTRONICS 

4.2.1  Outline  and  Block  Diagram 

When  a  bat  pulse  is  transferred  to  the  magnetic  drum  delay  unit  from  the  magnetic  tape  re¬ 
producer,  several  timing  and  gating  operations  take  place.  Figure  17  is  a  block  diagram  of  the 
arrangement  used  to  gate  a  bat  signal  onto  the  magnetic  drum.  On  the  second  channel  of  the  tape 
reproducer  a  pulse  marker  signal  has  been  placed  0.1  to  0.2  sec  before  the  start  of  the  bat  signal. 
This  signal  is  used  to  trigger  the  B  sweep  of  a  Tektronix  535A  oscilloscope.  The  oscilloscope 
is  used  as  a  signal  generator  to  provide  the  gating  signals  needed  to  gate  the  bat  signal  onto  the 
time  delay  drum.  The  A  sweep  on  the  oscilloscope  is  triggered  by  the  internal  delay  pickoff  on 
the  B  sweep.  The  A  sweep  gate  signal  is  used  to  drive  a  relay  that  connects  the  output  of  the 
record  amplifiers  of  the  Goff  delay  drum  to  the  record  heads.  The  A  sawtooth  signal  is  fed 
into  a  trigger  pickoff  unit  that  triggers  an  electronic  switch  on  and  then  off,  electronically  gating 
the  bat  signal  into  the  record  amplifier. 


SHIFT  CHANNEL  "B"  TRIGGER 

Fig.  17.  Electronics  used  to  record  signal  onto  magnetic  drum. 
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During  this  operation  of  the  record  electronics,  the  bat  signal  to  be  analyzed  is  played  by 
the  tape  reproducer.  The  gating  circuitry  gates  the  bat  signal  onto  track  1  of  the  drum.  The 
playing  speed  of  the  tape  reproducer  is  changed  by  the  incremental  amount  corresponding  to  the 
desired  amount  of  Doppler  shift;  then  the  tape  is  replayed  and  the  Doppler-shifted  bat  signal  is 
gated  onto  track  2  of  the  drum. 

4.2.2  Tape  Input 

The  slowed-down  bat-signal  sound  track  is  recorded  on  £-inch  magnetic  tape  at  l\  inches/sec 
on  a  conventional  stereo  tape  recorder,  an  Ampex  601-2.  After  locating  the  physical  position  of 
the  desired  bat  pulse  on  the  tape,  the  tape  is  magnetically  marked  as  follows: 

(a)  The  bat  pulse  is  manually  positioned  approximately  1  inch  ahead  of  the  record 
head,  with  the  tape  disengaged  from  the  capstan  drive  mechanism. 

(b)  The  recorder  control  is  turned  to  RECORD  position  and  then  off. 

This  procedure  magnetically  marks  both  tracks  of  the  magnetic  tape  at  the  points  where  the  re¬ 
cord  head  and  the  erase  head  are  in  contact  with  the  tape.  When  the  tape  is  played,  two  short 
pulses  are  reproduced  on  both  tracks  0.1  to  0.3  sec  before  the  bat  signal  is  reproduced  on  the  sig¬ 
nal  track  as  is  shown  in  Fig.  18. 


Fig.  18.  Synchronizing  marker  timing  on  typical 
1/4-inch  magnetic  tape  record  (bat  pulse  4). 
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The  first  marker  signal  on  the  control  track,  the  record  head  signal,  is  used  to  trigger  the 
gating  circuitry  that  gates  the  bat  signal  onto  the  magnetic  drum. 

Doppler  shift  is  introduced  into  the  tape-recorded  bat  signal  by  changing  the  speed  of  the  tape 
reproducer.  The  Ampex  601-2  tape  recorder  has  a  synchronous  motor  capstan  drive  so  that,  if 
the  frequency  of  the  motor  drive  power  is  changed,  the  speed  of  the  motor,  and  hence  the  tape 
speed,  will  be  changed  by  a  proportional  amount.  The  motor  of  the  Ampex  is  driven  by  a  variable 
frequency  oscillator  followed  by  a  high-power  amplifier.  A  very  stable  and  accurately  calibrated 
oscillator  is  needed,  since  the  change  in  frequency  due  to  the  Doppler  shift  is  usually  very  small. 
A  Krohn-Hite  model  310-AB  push-button  oscillator  was  found  to  be  suitable.  It  has  an  accuracy 
of  1  percent  and  a  drift  of  less  than  0.05  percent/hr.  The  power  amplifier  is  a  Vectron  VS-250 
that  was  modified  to  accept  the  external  oscillator  input. 

The  oscillator  setting  for  each  particular  Doppler  shift  was  calculated  by  using  the  approxi¬ 
mate  formula  for  the  oscillator  setting: 

f  =  6o  ( 1  —  -^r)  ■  (4-D 

A  list  of  the  oscillator  settings  for  various  target  velocities  is  given  in  Table  I. 
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TABLE  1 

DOPPLER  FACTOR  AND  MOTOR-DRIVE  OSCILLATOR  SETTINGS 
FOR  GIVEN  TARGET  VELOCITIES 

Target  Velocity 

Doppler  Factor 

Oscillator  Frequency 

(m/sec) 

(dimensionless) 

(cps) 

-10 

1.0579 

63.47 

—  8 

1.0463 

62. 78  \ 

-6 

1.0347 

62.09 

-4 

1.0232 

61.39 

-2 

1.0116 

60.70 

0 

1 . 0000 

60.00 

2 

0.9885 

59.31 

4 

0.  9769 

58.62 

6 

0.  9654 

57.  92 

8 

0.  9538 

57.  23 

10 

0.  9422 

56.54 

4.2.3  Electronic  Gating  and  Switching 

Two  types  of  gating  must  be  provided  to  record  a  chosen  bat  signal  onto  the  magnetic  drum. 
The  signal  is  electronically  switched  to  provide  the  desired  bat  pulse  to  the  record  amplifiers. 

The  currents  to  the  recording  heads  of  the  tape  drum  are  switched  so  that  recording  takes  place 
only  during  one  revolution  of  the  magnetic  drum.  The  electronic  switching  is  done  by  a  Grason- 
Stadler  (Model  E3382)  modulating  switch  and  the  head-current  switching  is  done  with  a  relay. 
Timing  signals  for  the  switching  processes  are  obtained  from  the  Tektronix  oscilloscope.  The 
signal  from  the  second  track  of  the  tape  reproducer,  after  being  amplified  and  filtered,  is  used 
to  trigger  the  B  sweep  of  the  Tektronix.  The  A  sweep  of  the  oscilloscope  is  triggered  internally 
after  a  fixed  time  delay.  The  timing  of  this  operation  is  shown  for  a  typical  bat  pulse  in  Fig.  19. 
The  A  sweep  gate  output  of  the  oscillator  is  used  to  drive  a  relay-driver  circuit  which,  in  turn, 
connects  the  record  amplifiers  to  the  record  heads  during  the  time  that  the  A  sweep  gate  is 

1  -45-7<9ll 


Fig.  19.  Timing  sequence  for  gating  signals  which  are  derived  from  oscilloscope 
that  is  triggered  by  channel  2  signal  of  input  tape. 
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positive.  The  A  sweep  sawtooth  is  used  to  drive  a  trigger  pickoff  circuit  which  turns  the 
electronic  switch  on  and  off  at  times  that  correspond  to  voltage  levels  on  the  sawtooth. 

The  relay-driver  circuit  provides  power  to  the  relay  when  the  oscilloscope  gate  is  positive, 
that  is  when  the  gate  is  at  +20  volts.  This  gate  is  usually  set  for  a  duration  of  between  0.2  to 
0.25  sec  because  the  drum  makes  one  revolution  in  0.25  sec.  The  schematic  for  the  relay  driver 
is  shown  in  Fig.  20.  The  relay  ungrounds  the  output  of  the  record  amplifiers  and  connects  the 
amplifiers'  outputs  to  their  respective  record  heads  or  to  ground  depending  on  whether  the  record 
switch  is  on  or  off.  In  Goff' s  design  of  the  magnetic  drum  circuitry,  the  movable  head  is  a  record 
head.  After  the  pulse  signals  have  been  recorded  on  the  drum,  it  is  one  of  the  playback  heads 
that  must  be  moved.  Fortunately,  the  movable  head  works  satisfactorily  as  a  playback  head  and 
its  connections  have  therefore  been  modified  so  that  it  can  be  used  as  the  playback  head  as  well 
as  the  record  head.  A  switch  has  been  installed  on  the  input  of  the  number  2  playback  amplifier 
to  facilitate  the  use  of  either  the  normal  record  or  the  normal  playback  head. 

The  high-frequency  bias  signal  to  the  erase  heads  is  also  switched.  Before  recording  a  pulse 
on  the  drum,  the  entire  drum  is  erased  by  turning  on  the  erase  switch.  If  this  switch  is  turned 
off,  a  large  spike  will  be  recorded  on  the  magnetic  drum.  To  avoid  this,  the  bias  to  the  erase 
heads  is  turned  off  slowly  by  disconnecting  B+  from  the  bias  oscillator.  This  is  accomplished  by 
turning  off  the  microswitch  that  provides  power  to  a  relay  in  the  record  amplifier  unit  which,  in 
turn,  switches  the  power  to  the  drum  drive  motor  as  well  as  switching  the  B+  to  the  bias  oscillator. 
After  turning  off  the  bias  in  this  way,  the  erase  switch  is  turned  off  and  then  the  microswitch  is 
turned  back  on.  The  entire  drum  is  now  erased. 

The  A  sweep  sawtooth  from  the  oscilloscope  is  a  positive-going  ramp  starting  at  ground 
potential  and  climbing  to  about  +150  volts.  Pulses  to  turn  the  electronic  switch  on  and  off  are  ob¬ 
tained  from  a  trigger  pickoff  circuit  which  operates  from  the  A  sweep  sawtooth.  The  schematic 
diagram  for  the  trigger  pickoff  circuit  is  given  in  Fig.  21.  The  50-volt  negative  pulses  needed  to 
switch  the  electronic  switch  are  obtained  from  the  sawtooth  with  Philbrick  operational  amplifiers 
(type  K-2W).  The  transition  time  on  the  output  of  these  amplifiers  is  decreased  by  the  use  of 
positive  feedback  around  the  amplifiers.  The  triggering  levels  for  turnon  and  turnoff  are  set  by 
the  potentiometers  on  the  positive  input  of  the  operational  amplifiers.  Figure  19  shows  the  timing 
of  the  turnon  and  turnoff  signals  for  a  typical  bat  pulse.  When  the  voltage  on  the  negative  input  to 
the  amplifier  increases  to  equal  the  voltage  on  the  positive  input,  the  output  goes  through  a  fast 
negative  transition.  This  negative  step  is  capacitor-coupled  to  the  triggering  input  of  the  Grason- 
Stadler  switch.  The  diode  in  parallel  with  the  output  is  necessary  because  of  the  peculiarities  of 
the  Grason-Stadler  switch.  If  the  trigger  signal  to  this  switch  is  allowed  to  go  negative,  the  unit 
will  not  turn  off  properly. 

4.2.4  Drum  Delay  Unit 

The  drum  delay  unit  was  originally  designed  to  operate  by  continuously  recording  and  re¬ 
producing  signals  on  two  separate  tracks.  It  has  been  modified  to  accept  pulse  signals  that  are 
recorded  once  and  reproduced  many  times. 

The  drum  is  8  inches  in  diameter  and  rotates  4  revolutions/sec,  giving  a  linear  speed  of 
100  inches/sec  at  the  rim  of  the  drum.  The  record  and  reproduce  heads  are  located  a  nominal 
0.001  inch  from  the  surface  of  the  drum.  For  an  input  signal  of  constant  amplitude,  the  recording 
current  is  constant  to  about  1  kcps  and  gradually  increases  to  a  maximum  of  pre-emphasis  of  7  db 
at  10  kcps.  The  reproduce  amplifiers  have  several  filters  that  can  be  switched  to  give  the  frequency 
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Fig.  20.  Relay  driver  and  head  switching. 
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characteristics  that  are  shown  in  Fig.  22.  For  the  work  with  the  bat  signals,  position  three  was 
normally  used  on  the  reproduce  amplifiers  to  help  eliminate  the  small  amount  of  60-cps  hum  that 
is  introduced  by  the  drum  delay  unit. 

The  drum  surface  is  magnetically  very  smooth  and  there  do  not  seem  to  be  any  dropouts  in 
signals  that  are  recorded  on  the  drum.  Because  of  a  very  slight  eccentricity  in  the  bearings  that 
hold  the  movable  head,  the  drum  output  varies  as  a  function  of  the  position  of  the  movable  head. 
The  gains  of  the  drum  amplifiers  are  not  stabilized  so  a  slow  variation  in  the  output  of  the  drum 
is  also  encountered.  The  addition  of  attenuators  on  the  output  of  the  reproduce  electronics  helps 
to  compensate  for  these  variations.  After  the  pulse  signals  are  recorded  on  the  drum,  these  at¬ 
tenuators  are  set  so  that  the  output  is  within  0.2db  of  the  nominal  output,  3  volts  peak  to  peak. 

A  small  amount  of  noise  is  introduced  into  the  system  by  a  permanently  recorded  signal  on 
the  drum.  This  signal  is  probably  due  to  variations  in  the  magnetization  of  the  magnetic  layer 
in  the  vicinity  of  the  recorded  track  that  is  too  far  out  to  be  erased  by  the  erase  head.  Goff  men¬ 
tions  in  his  instruction  manual  (see  Ref.  17,  p.  41)  that  such  a  problem  exists  with  this  equipment 
and  says  that  this  signal  can  be  erased  with  a  hand-held  erase  head,  an  oscillator,  and  a  power 
amplifier.  Since  this  signal  is  about  40  db  below  the  level  of  the  recorded  signal,  the  process  of 
hand  erasing  was  not  deemed  necessary. 
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Fig.  21.  Trigger  pickoff  circuit. 
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Fig.  22.  Goff  time  delay  drum:  frequency  response  for  various  positions  of  low-frequency 
rolloff  switch.  (Reproduced  from  Goff's  "Instruction  Manual  for  Analog  Correlator,"  1954, 
with  permission  from  John  A.  Kessler,  former  Executive  Officer  of  the  M.  I.  T.  Acoustics 
Laboratory.) 

The  drum  unit  has  a  mechanical  filter  to  help  reduce  flutter  in  the  reproduced  signal.  The 
flutter  of  the  unit  is  below  0.02  percent.  The  drum  is  driven  by  a  synchronous  motor  that  is,  in 
turn,  driven  by  a  precision-fork  frequency  standard  through  a  power  amplifier. 

The  movable  head  can  be  moved  manually  or  driven  at  various  rates  by  a  synchronous  motor- 
and-gear  arrangement.  The  head  has  access  to  205  msec  of  the  250-msec-long  surface  of  the 
drum.  There  is  a  mechanical  revolution  counter  on  a  shaft  that  is  geared  to  the  head  drive  shaft. 
This  counter  displays  the  drum  delay  in  tenths  of  a  millisecond.  Since  the  signals  are  put  on  at 
random,  this  counter  can  indicate  only  a  relative  position  of  the  movable  head. 

A  potentiometer  with  a  slip-clutch  drive  was  added  to  the  drum  delay  unit  to  give  an  electri¬ 
cal  output  that  is  proportional  to  the  position  of  the  movable  head.  Figure  23  shows  how  the  ten- 
turn  potentiometer  is  mounted  on  the  drum.  The  potentiometer  shaft  is  gear-driven  from  the 
motor  that  drives  the  movable  head.  A  pair  of  gears  can  be  changed  to  give  various  ratios  be¬ 
tween  the  drum  delay  time  and  the  potentiometer  rotation,  as  shown  in  Table  II.  The  electrical 
circuitry  and  the  performance  characteristics  of  the  potentiometer  position  indication  are  dis¬ 
cussed  in  Sec.  4.3.4. 


Fig.  23.  Mounting  of  ten-turn  heliopot  that 
is  used  to  drive  X-axis  of  X-Y  recorder. 
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TABLE  II 

DRUM  DELAY  TIME/ROTENTIOMETER-ROTATION  CALIBRATION 

Gear  Ratio 

Calibration 

Drum  Delay  Time/Potentiometer  Rotation 

Motor  Gear/Potentiometer  Gear 

(msec/ revo  1  ut  i  on) 

1:4 

0.625 

1:2 

1.25 

1:1 

2.5 

2:1 

5.0 

4:1 

10.0 

4.3  REPRODUCE  ELECTRONICS 
4.3.1  Outline  and  Block  Diagram 

Once  the  signal  is  on  the  magnetic  drum,  the  problem  becomes  that  of  multiplying  and  inte¬ 
grating  the  function  given  by  Eq.  (2-27): 


x(t,  a) 


I 


u(t)  u[a(t  —  t)]  4a  dt 


[  (2-27)] 


After  taking  care  of  the  time  delay  and  the  Doppler  shift,  the  only  operations  left  to  perform  are 
the  multiplication  and  the  integration.  Multiplication  by  4~a  is  omitted,  since  the  error  intro¬ 
duced  by  omitting  this  multiplication  is  less  than  3  percent.  A  block  diagram  of  the  reproduce 
electronics  is  given  in  Fig.  24.  The  multiplication  operation  is  performed  with  a  tube- 
characteristic  quarter-square  multiplier.  The  product  signal  is  gated  and  then  integrated  by 
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Fig.  24. 


nics  used  for  multiplication  and  integration  of  signals  on  magnetic  drum. 
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using  relays  and  Philbrick  operational  amplifiers.  The  output  of  the  integrator  is  sampled  at  the 
end  of  the  integration  and  is  electrically  held  on  a  capacitor  in  a  sample-and-hold  circuit.  The 
output  of  the  sample -and-hold  circuit  is  connected  to  an  amplifier  that  drives  the  Y  input  of  an 
X-Y  recorder. 

The  X  drive  of  the  X-Y  recorder  is  driven  from  either  a  12-minute  sawtooth  generator  or 
from  the  output  of  an  amplifier  connected  to  the  potentiometer  that  is  mounted  on  the  drum  delay 
unit.  Since  the  temporal  stability  and  the  linearity  of  the  potentiometer  drive  are  much  better 
than  that  of  the  sawtooth  generator,  the  sawtooth  generator  is  no  longer  used. 

When  three-dimensional  drawings  are  to  be  made  from  the  output  of  the  computer,  it  is  de¬ 
sirable  to  plot  the  slices  of  the  ambiguity  diagram  on  a  30°  angle  with  respect  to  the  horizontal. 
This  is  accomplished  by  algebraically  subtracting  tan  30°  times  the  X  input  to  the  X-Y  recorder 
from  the  output  of  the  sample-and-hold  circuit  and  applying  this  to  the  Y  input  of  the  X-Y  recorder. 

4.3.2  Multiplier 

The  multiplier  is  a  direct-coupled  multiplier  that  uses  the  nearly  parabolic  plate-current- 
grid-voltage  characteristics  of  back-to-back  6AS6/57  25  pentodes  with  the  control  and  suppressor 
grids  tied  together  in  a  quarter-square  arrangement.  The  output  of  the  multiplier  is  ideally 
given  by 

xy  =  j  ((x  +  y)2  -  (x  -  y)2]  .  (4-2) 

The  multiplier  consists  of  a  pair  of  difference  amplifiers  on  the  two  inputs,  a  resistive  summing 

network,  four  6AS6/5725  squaring  tubes,  and  a  difference  amplifier  on  the  output.  All  amplifiers 

are  direct  coupled,  giving  the  multiplier  a  DC  to  better  than  lOOkcps  frequency  response.  Since 
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the  multiplier  is  written  up  in  some  detail  in  the  Master's  thesis  report,  only  the  changes  and 
the  shortcomings  will  be  mentioned  here. 

A  very  serious  amount  of  DC  drift  was  encountered  in  the  output  of  the  original  multiplier 
Therefore,  a  new  power  supply  was  designed  and  built  into  the  multiplier  which  made  all  the  power 
supply  voltages  dependent  on  one  accurate  gas-tube  standard.  This  modification  decreased  the 
DC  drift  by  about  a  factor  of  two,  not  nearly  enough.  It  was  found  that  a  variation  in  the  voltage 
of  the  power  line  caused  a  sizeable  change  in  the  output  level.  This  effect  is  probably  due  to  the 
fact  that  the  heaters  in  the  tubes  are  not  run  from  a  regulated  supply.  With  the  addition  of  a  Sola 
constant -voltage  transformer  on  the  power  input  of  the  multiplier,  thus  regulating  the  heaters, 
the  output  drift  due  to  changes  in  the  line  voltage  decreased  by  more  than  an  order  of  magnitude. 

Normally  the  multiplier  operates  with  1-volt  rms  signals  on  the  input.  The  output  for  this 
input,  with  different  frequency  sine  waves  on  the  two  inputs,  is  about  2  volts  rms.  After  the 
above  changes  were  made  to  the  multiplier,  the  output  drift  was  on  the  order  of  0.01  volt. 

The  resistors  in  the  voltage  dividers  that  feed  the  output  difference  amplifier  change  value 
slightly  with  a  change  in  temperature.  Blowing  cold  air  on  one  of  these  resistors  changes  the 
output  level  by  about  0  1  volt.  Wire-wound,  deposited  film,  and  common  carbon  resistors  were 
tried  in  this  location  without  causing  a  significant  change  in  the  amount  of  drift.  Since  the  re¬ 
sistors  are  one-watt  resistors  with  only  l/8  watt  of  power  dissipated,  it  did  not  seem  that  an 
increase  in  the  wattage  would  give  a  substantial  decrease  in  the  amount  of  drift  encountered. 

Since  the  0.01 -volt  drift  is  still  greater  than  the  resolution  capabilities  of  the  output  re¬ 
corder,  a  circuit  for  effectively  resetting  the  output  level  of  the  multiplier,  before  multiplication 
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of  each  pulse,  was  devised.  This  system  ties  the  two  inputs  of  the  integrator  circuit  together 
before  the  pulse  comes  through  the  multiplier.  The  integral  is  then  the  integral  of  the  difference 
between  the  multiplier  outputs  with  the  pulse  on  the  inputs  and  with  zero  on  the  inputs.  This  sys¬ 
tem  satisfactorily  solves  the  DC  drift  problem  and  is  described  in  more  detail  in  Sec.  4.3.3. 

4.3.3  Pulse  Integration  and  Control  Circuits 

Gating,  integration,  and  sampling  are  done  with  Philbrick  K-2W  operational  amplifiers  con¬ 
trolled  by  relays.  The  problem  is  one  of  integrating  a  portion  of  the  signal  from  the  multiplier, 
holding  the  value  of  the  integral  for  display,  and  resetting  the  integrator  to  prepare  for  the  next 
pulse  from  the  multiplier. 

The  product  signal  from  the  multiplier  is  amplified  by  a  factor  of  ten  before  gating  opera¬ 
tions  take  place  so  that  its  amplitude  is  in  the  optimum  range  for  the  operational  amplifiers. 

This  is  done  with  suitable  feedback  to  give  the  gain  of  minus  ten.  Before  the  arrival  of  the  prod¬ 
uct  pulse,  the  plus  and  minus  inputs  to  the  integrator  are  tied  together.  Just  before  the  pulse 
arrives,  the  positive  input  to  the  integrator  is  disconnected  from  the  output  of  the  amplifier. 

There  is  a  1-fif  capacitor  holding  the  potential  on  the  positive  input  of  the  integrator.  When  the 
pulse  signal  arrives,  the  difference  between  the  signal  and  the  potential  stored  on  the  capacitor 
is  integrated.  Since  the  capacitor  was  set  when  the  multiplier  was  multiplying  zero  by  zero,  the 
output  of  the  integrator  is  the  integral  of  the  difference  between  the  product  signal  from  the  mul¬ 
tiplier  and  the  zero  level  of  the  product  signal. 

After  integration,  the  output  of  the  integrator  is  sampled  by  the  sample-and-hold  circuit. 
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This  sample-and-hold  circuit  is  similar  to  the  one  described  by  Brubaker.  It  consists  of  an 
operational  amplifier  with  feedback  for  a  gain  of  plus  one  followed  by  an  integrator,  connected 
through  a  sampling  switch,  with  over-all  feedback  on  the  circuit  to  give  a  gain  of  minus  one  when 
the  sampling  switch  is  closed.  At  the  end  of  the  integration  by  the  integrator  circuit,  the  sam¬ 
pling  relay  is  closed  for  about  10  msec.  The  contacts  on  the  sampling  relay  bounce  for  about 
2  msec  after  they  first  close.  The  output  of  the  sampling  circuit  rises  (or  falls)  in  a  linear  fash¬ 
ion,  while  the  contacts  are  closed,  until  it  is  within  about  2  volts  of  the  negative  of  the  input  volt¬ 
age.  The  output  then  exponentially  approaches  the  negative  of  the  input  voltage  with  approximately 
a  0.25-msec  time  constant. 

When  the  echolocation  pulse  is  recorded  on  the  magnetic  drum,  the  record  head  current  is 
switched  on  and  then  off.  This  switching  of  the  head  current  causes  two  large  pulses  to  be  re¬ 
corded  on  the  magnetic  drum.  One  of  these  pulses  is  used  to  trigger  a  string  of  one-shot  multi¬ 
vibrators  that  control  the  timing  of  relay  operations  in  the  pulse  integrator.  The  signal  from 
channel  one  of  the  magnetic  drum  is  fed  into  a  Schmidt  trigger  circuit  that  shapes  this  trigger 
pulse.  This  shaped  trigger  pulse  sets  a  flip-flop  which,  in  turn,  triggers  the  delay  one-shot 
multivibrator.  The  timing  of  the  operations  of  this  circuit  is  shown  in  Fig.  25.  When  the  delay 
one-shot  turns  off,  the  integrate  one-shot  is  triggered,  operating  the  integrate  relay  which  dis¬ 
connects  the  positive  input  of  the  integrator  from  the  output  of  the  gain-of-ten  amplifier.  The 
integrator  circuit  is  now  integrating  the  output  of  the  multiplier.  When  the  integrate  one-shot 
turns  off,  the  sample  one-shot  is  triggered,  the  sample  relay  closes,  and  the  integrate  relay 
opens.  The  capacitor  in  the  sample-and-hold  circuit  is  charged  to  the  proper  value  and  then, 
after  10  msec,  the  sample  one-shot  turns  off.  The  reset  one-shot  is  triggered  and  the  sample 
relay  drops  out;  then  the  reset  relay  closes,  discharging  the  integration  capacitor  in  the  inte¬ 
grator,  and  remains  closed  until  the  reset  one-shot  turns  off.  When  the  reset  one-shot  turns  on. 
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Fig.  25.  Timing  of  the  sequence  of  operations  performed 
by  the  pulse  integrator  timing  circuitry. 

the  flip-flop  is  reset.  Until  this  flip-flop  is  reset  and  the  reset  one-shot  turns  off,  the  delay 
one-shot  cannot  be  triggered.  This  prevents  the  echolocation  signal  from  triggering  the  string 
of  one-shot  multivibrators. 

Therefore,  one  trigger  signal  is  propagated  down  the  string  of  one-shot  multivibrators  and 
then  returns  to  enable  the  flip-flop.  If,  while  setting  the  delays  of  the  one-shot  multivibrators, 
the  trigger  signal  is  lost  (this  can  and  does  happen  under  certain  circumstances  when  the  one- 
shot  capacitors  are  being  switched),  the  flip-flop  can  be  reset  manually  by  pressing  the  reset 
push  button. 

4.3.4  X-Y  Plotter  Output  and  Associated  Electronics 

A  Variplotter,  Model  1100E,  X-Y  pen-and-ink  recorder  is  used  to  plot  the  output  of  the 
pulsed  integrator  as  a  function  of  the  delay  time.  Since  the  movable  head  on  the  delay  drum  can 
be  driven  at  a  constant  rate  by  a  synchronous  motor,  it  was  decided  that  the  easiest  way  to  drive 
the  X  input  of  the  X-Y  recorder  was  from  a  sawtooth  generator  that  could  be  started  at  the  same 
time  that  the  synchronous  motor  was  started.  A  12-minute  sawtooth  generator  was  designed  and 
constructed.  This  generator  uses  a  100-megohm  resistor  and  five  1-pf  capacitors  in  an  inte¬ 
grator  circuit.  The  integrator  is  a  standard  integrator  built  around  a  Philbrick  K-2W  opera¬ 
tional  amplifier.  A  relay  is  used  to  reset  the  sawtooth  generator  by  shorting  out  the  integrating 
capacitor.  The  output  sawtooth  sweeps  from  —40  volts  up  to  +35  volts.  This  is  enough  to  give 
the  required  15  inches  of  deflection  on  the  recorder  when  the  recorder  sensitivity  is  set  at 
5  volts/inch. 

After  a  considerable  amount  of  computer  operation,  it  was  decided  that  the  12-minute  saw¬ 
tooth  generator  could  stand  some  improvement.  Although  this  generator  was  accurate  to  better 
than  one  percent,  it  had  the  wrong  type  of  accuracy.  Since  a  plot  is  made  from  a  minus  time 
delay  through  zero  up  to  a  plus  time  delay,  even  though  the  minus  time  delay  is  set  exactly,  the 
position  of  the  zero  might  vary  back  and  forth  1  percent  of  the  full-scale  time  delay.  If  errors 
are  to  be  introduced,  it  is  much  better  that  they  are  introduced  at  the  two  extremes  of  the  time 
delay  rather  than  have  the  zero  displaced. 
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A  method  of  marking  the  zero  of  the  graph  was  used.  When  a  monitor  oscilloscope  registered 
that  the  zero  point  had  been  reached,  a  button  was  pressed,  causing  a  l/Z -  to  1-inch  jump  in  the 
pen  position.  The  accuracy  of  pressing  the  button  at  the  point  of  zero  time  delay  was  approxi¬ 
mately  l/4  to  l/2  sec,  corresponding  to  0.005  to  0.010  inch  on  the  output  graph.  The  only  trouble 
encountered  in  pressing  the  button  was  caused  by  the  fact  that  the  button  was  pressed  6  minutes 
after  the  start  of  the  plotting,  halfway  through  the  12-minute  plot.  Operator  fatigue  caused  the 
marking  point  to  be  missed  several  times.  A  photographic  timer  was  used  very  successfully  in 
alerting  the  operator.  This  eliminated  the  marking  being  missed  due  to  operator  fatigue. 

Another  solution  of  the  problem  of  marking  the  zero  point  is  to  use  a  potentiometer  to  obtain 
the  X  deflection  of  the  X-Y  recorder. 

The  potentiometer  is  mounted  on  the  drum  unit,  gear  driven  from  a  shaft  that  is  driven  by 
the  synchronous  motor  that  drives  the  movable  head.  A  mechanical  description  of  the  mounting 
and  gearing  of  this  potentiometer  is  given  in  Sec.  4.  2.  4.  The  resistors  in  parallel  with  the  po¬ 
tentiometer  that  go  to  ground  help  to  reduce  the  effect  of  possible  power  supply  variations  on  the 
nulled  output  of  the  potentiometer.  This  circuit  keeps  the  zero  of  the  delay  scale  of  the  output 
within  an  ink-line  width,  0.02  inch,  of  the  correct  position,  which  is  more  than  an  order  of  mag¬ 
nitude  better  than  the  sawtooth  generator.  Not  only  does  the  potentiometer  drive  improve  the 
accuracy  of  the  location  of  the  zero  position  of  the  graph,  but  it  also  improves  the  linearity  of 
the  sweep  and  the  temporal  stability  of  the  delay  scale  calibration.  The  linearity  of  the  sawtooth 
generator  and  of  the  potentiometer  circuit  is  illustrated  in  Fig.  26(a-b).  The  error  from  linearity 
is  given  as  a  function  of  the  desired  X  position  of  the  pen  on  the  X-Y  recorder.  This  error  from 
linearity  is  plotted  for  four  consecutive  sweeps  by  each  sweep  system.  It  may  be  readily  observed 
that  use  of  the  potentiometer  circuit  has  greatly  improved  the  temporal  stability. 

When  a  three-dimensional  drawing  is  to  be  made  from  the  computer  output,  it  is  desirable 
to  plot  the  slices  of  the  ambiguity  surface  on  a  slant.  A  typical  slice  is  shown  in  Fig.  27.  A  tri¬ 
metric  representation  was  chosen  where  one  axis  is  plotted  on  a  30°  slope  with  respect  to  the 
horizontal  and  the  other  axis  is  plotted  on  a  10°  slope  with  respect  to  the  horizontal.  These  slices 
are  easily  obtained  by  simply  subtracting  tan  30°  times  the  X  deflection  of  the  recorder  from  the 
Y  deflection  of  the  recorder.  The  circuit  that  does  this  is  included  on  the  X-Y  recorder-drive 
chassis.  The  signal  from  the  sample-and-hold  circuit  is  the  negative  of  the  cross-correlation 
function,  so  that  when  the  X  deflection  signal  is  added  to  it  at  the  summing  point  of  the  opera¬ 
tional  amplifier  with  the  appropriate  multiplicative  constant,  the  output  of  the  operational  am¬ 
plifier  is  equal  to  the  cross-correlation  function  minus  tan  30°  times  the  X  deflection.  An  atten¬ 
uator  is  inserted  between  the  sample-and-hold  circuit  and  the  summing  operational  amplifier  to 
adjust  the  gain  of  the  system  so  that  the  maximum  amplitude  of  the  autocorrelation  function  can 
be  normalized  to  one. 

4.4  ERRORS  AND  LIMITATIONS  IN  THE  ACCURACY  OF  THE 
AMBIGUITY  DIAGRAM  COMPUTER 

This  section  pertains  to  an  investigation  of  the  errors  introduced  by  the  ambiguity  diagram 
computer.  Each  significant  error  is  discussed  in  an  attempt  to  place  some  meaningful  bounds 
on  the  accuracy  of  the  computer  output. 

4.4.1  Function  Error 

An  error  in  the  functions  u(t)  and  u(a(t  —  t))  can  be  introduced  in  several  places.  The  proc¬ 
ess  of  recording  the  bat  signal  onto  tape  is  probably  the  place  where  the  largest  amount  of  error 


43 


POSITION  ERROR  (in.)  POSITION  ERROR  (in.) 


1 3-45-T496(o -b)l 


(a)  12-minute  sawtooth  generator. 


(b)  Potentiometer  drive. 


Fig.  26.  Linearity  of  sweep  output. 


Fig.  27.  Slice  of  typical  ambiguity  surface  plotted  on  30° slope  for  use  in  making  trimetric 
representation  of  surface  (bat  pulse  4,  a  =  1). 
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is  introduced.  The  microphones  that  are  used  to  pick  up  the  bat  signal  are  not  perfect  and  have 
a  frequency-dependent  gain  pattern.  Echoes  from  the  surroundings  and  noise  from  the  micro¬ 
phone  amplifiers  add  a  level  of  hash  to  the  signals.  It  is  difficult  to  estimate  in  some  cases  which 
is  the  white  noise  introduced  by  the  microphone  amplifiers  and  which  is  the  noise  introduced  by 
room  echoes.  This  hash  level  can  be  as  low  as  —  35  db,  and  is  usually  not  higher  than  —  20  db. 

The  signal  that  is  recorded  on  the  magnetic  tape  will  therefore  be  assumed  to  be  the  signal  for 
which  the  ambiguity  diagram  is  to  be  calculated. 

When  the  tape  is  played,  a  variation  in  the  amplitude  of  the  signal  might  occur.  With  the 
Ampex  601-2  tape  recorder,  using  type  ^ -200-24  Minnesota  Mining  and  Manufacturing  magnetic 
tape,  the  amplitude  variation  appeared  to  be  less  than  1  percent. 

When  the  signal  is  put  onto  the  magnetic  drum,  the  drum  noise  is  added  to  the  signal.  The 
level  of  the  noise  on  the  magnetic  drum  is  about  40  db  below  the  level  of  the  recorded  signal. 

This  corresponds  to  a  signal  error  of  1  percent. 

4.4.2  Multiplier  Error 

The  error  introduced  by  the  multiplier  is  probably  the  most  difficult  to  estimate.  There  can 
be  several  types  of  errors  introduced  by  the  multiplier.  Most  of  these  do  not  affect  the  results 
of  the  correlation  integral.  If,  for  example,  there  is  a  drift  in  the  DC  balance  of  the  input  dif¬ 
ference  amplifiers,  the  output  of  the  mulitplier  would  be  given  by 

eQ  =  (e1  +  a)  (e2  +  b)  (4-3) 

where  e^  and  e^  are  the  input  voltages  to  the  multiplier  and  a  and  b  are  the  unbalance  levels. 
Equation  (4-3)  can  be  rewritten  as 

eQ  =  e^e^  +  be1  +  ae^  +  ab  (4-4) 

The  last  term  of  this  equation  does  not  affect  the  output  since  this  corresponds  to  a  steady  DC 
level  on  the  output  of  the  multiplier.  This  steady  level  is  canceled  out  by  the  pulsed  integrator 
circuit  described  in  Sec.  4.3.3.  The  two  middle  terms  in  Eq.  (4-4)  do  not  affect  the  output  of  the 
integrator,  since  neither  the  e^  nor  the  e^  signal  has  any  DC  component  in  its  energy  density 
spectrum. 

With  two  sine  waves  on  the  input  to  an  ideal  multiplier,  the  output  contains  only  the  sum  and 
difference  frequencies.  With  the  multiplier  used  for  these  experiments,  the  output  contains  not 
only  the  sum  and  difference  frequencies,  but  also  components  of  the  two  input  frequencies  and 
cross-modulation  terms  between  the  output  and  the  sum  frequencies.  The  level  of  these  spurious 
modulation  products  was  usually  at  least  20  db  below  the  level  of  the  sum  and  difference  frequency 
components. 

Since  the  only  term  that  appears  on  the  output  of  the  integrator  is  the  zero  frequency  compo¬ 
nent,  the  testing  of  the  multiplier  accuracy  is  done  only  on  the  difference  component  of  the  output. 
A  block  diagram  of  the  testing  arrangement  is  shown  in  Fig.  28.  Two  oscillators  feed  the  two  in¬ 
puts  of  the  multiplier  through  attenuators.  The  output  of  the  multiplier  is  low -passed  and  its  level 
is  measured  with  a  true-rms  voltmeter.  The  reading  on  the  voltmeter  is  plotted  as  a  function  of 
the  signal  level  on  the  X  input  to  the  multiplier  for  three  different  Y  input  levels  and  is  shown 
in  Fig.  29.  The  maximum  deviation  from  the  correct  output  for  a  50-db  operating  range  is  less 
than  0.5  db.  This  corresponds  to  an  error  of  less  than  0.4  percent. 
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Fig.  28.  Setup  used  to  test  accuracy  of  multiplier. 


Fig.  29.  Mu Itiplier  accuracy  test.  Three  straight  lines  indicate  correct 
output  for  multiplier;  circles  show  actual  output  of  multiplier. 
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INDICATED  FREQUENCY  (cps) 

Fig.  30.  Error  in  frequency  of  Doppler  shift  oscillator,  normalized 
relative  to  frequency  at  setting  of  60  cps. 
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4.4.3  Doppler  Error 


The  error  in  the  accuracy  of  the  Doppler  shift  puts  a  serious  limitation  on  the  accuracy  of 
the  ambiguity  diagrams.  The  largest  Doppler  shift  that  was  plotted  in  this  project  was  10  m/sec. 
This  corresponds  to  a  Doppler  factor  of  1.06,  or  in  other  words  is  only  a  6-percent  change  in  the 
frequency  of  the  bat  signal.  Additional  changes  in  the  frequency  of  the  bat  signal  can  be  intro¬ 
duced  in  the  following  three  ways. 

(a)  Oscillator  Error:—  The  Doppler  shift  oscillator  that  drives  the  tape  repro¬ 
ducer  can  introduce  errors  in  the  Doppler  shift  of  the  signal,  if  it  is  not 
operating  at  the  correct  frequency.  The  specifications  for  the  oscillator 
given  by  the  manufacturer  state  that  the  oscillator  has  an  accuracy  of  1  per¬ 
cent  and  a  drift  of  less  than  0.05  percent/hr.  The  absolute  accuracy  of  the 
oscillator  is  of  no  real  importance  for  this  work,  since  it  is  the  relative 
accuracy  of  the  oscillator  that  affects  the  accuracy  of  the  Doppler  shift. 

That  is  to  say,  if  the  oscillator  were  operating  at  a  frequency  1  percent 
above  the  desired  frequency  for  all  the  oscillator  settings,  no  change  in 

the  output  of  the  computer  would  result.  If,  however,  the  oscillator  oper¬ 
ates  at  1  percent  above  the  desired  frequency  for  one  setting  and  at 
0.96  percent  above  the  desired  frequency  for  another  setting,  a  0.04-percent 
error  in  the  position  of  the  Doppler  slice  would  result. 

Figure  30  illustrates  the  relative  amount  of  error  in  the  oscillator  fre¬ 
quency  for  the  various  settings  of  the  oscillator.  This  graph  has  been  nor¬ 
malized  in  such  a  way  that  the  error  in  the  oscillator  frequency  is  meas¬ 
ured  relative  to  the  oscillator  frequency  that  occurs  when  the  oscillator  is 
set  on  60cps.  It  is  obvious  from  this  graph  that  the  relative  error  in  the 
position  of  the  Doppler  slices  due  to  the  oscillator  is  less  than  0.05  percent. 

(b)  Tape  Reproducer  Flutter:—  The  tape  reproducer  can  introduce  errors  of 
two  types.  It  can  operate  at  a  speed  lower  or  higher  than  the  correct 
speed  or  it  can  actually  change  speed  during  the  time  that  the  bat  signal 

is  being  played.  The  manufacturer  specifies  that  the  flutter  (error  around 
the  operating  speed)  is  less  than  0.17  percent.  The  tape  reproducer  was 
analyzed  by  the  method  described  by  Cox^O  and  was  found  to  have  a  peak 
flutter  of  0.36  percent.  If  flutter  components  above  10  cps  are  ignored, 
the  flutter  is  less  than  the  0.17  percent  that  is  specified  by  the  manufacturer. 

The  time  displacement  error  accumulated,  in  the  l/4sec  that  it  takes  the 
drum  to  make  one  revolution,  is  0.2  msec.  Along  with  the  flutter,  the  tape  re¬ 
producer  introduces  an  operating  speed  error.  This  error  is  0.075  percent. 

(c)  Drum  Flutter:—  According  to  Goff  (see  Ref.  15,  p.  19),  the  amount  of  flutter 
introduced  by  the  magnetic  drum  is  less  than  0.02  percent.  This  has  been 
measured  and  confirmed. 

The  total  error  in  the  position  of  a  Doppler  slice  is  the  sum  of  the  individual  frequency 
errors  divided  by  the  percentage  change  in  frequency  that  corresponds  to  the  largest  Doppler 
shift.  Another  way  of  calculating  this  is  to  determine,  from  the  frequency  errors,  the  actual 
amount  of  velocity  that  can  be  expected.  The  maximum  frequency  error  due  to  the  flutter  of  the 
tape  reproducer  and  the  magnetic  drum  is  0.38  percent.  The  maximum  frequency  error  due  to 
oscillator  drift  and  the  operating  speed  error  of  the  tape  reproducer  is  0.125  percent.  Adding 
these  two  errors,  the  maximum  frequency  error  is  0.5  percent.  If  the  velocity  of  sound  in  air 
is  taken  as  346  m/sec,  this  frequency  error  corresponds  to  a  target  velocity  error  of  0.9  m/sec. 
If  the  ambiguity  diagram  extends  from  a  target  velocity  of  —10  m/sec  to  a  target  velocity  of 
+  10m/sec,  the  0.9  m/sec  is  a  position  error  of  the  Doppler  slice  equal  to  4.4  percent. 

4.4.4  Time  Delay  Error 

The  time  delay  error  caused  by  the  potentiometer  drive  system  and  the  X-Y  plotter  can  be 
estimated  from  Fig.  26(b).  Since  the  maximum  error  indicated  is  0.02  inch  and  the  full-scale 
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TABLE  III 

SUMMARY  OF  THE  ERRORS  INTRODUCED  BY  THE 
DIAGRAM  COMPUTER 

AMBIGUITY 

Error 

Change 

(percent) 

Total  Error 
of  Full  Scale 
(percent) 

Height 

Function  error 

4.0 

Tape  dropout  (channel  1) 

1.0 

Tape  dropout  (channel  2) 

1.0 

Drum  noise  (channel  1) 

1.0 

Drum  noise  (channel  2) 

1.0 

Multiplier 

0.4 

Integrator 

1.0 

Sample  and  hold 

0.00025 

5.4 

Time  delay 

Potentiometer  drive 

0. 14 

Goff  drum 

0.2 

Zero  setting 

0.2 

0.54 

Doppler  shift 

Osci  1  lator 

0. 05 

0.43 

Tape  reproducer 

Flutter 

0. 36 

3.  1 

Drift 

0. 075 

0.65 

Drum  flutter 

0.02 

0.  17 

4.4 
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width  of  the  X-Y  plotter  is  15  inches,  the  maximum  error  in  the  time  delay  is  estimated  at 
0.14  percent. 

The  error  in  the  calibration  of  the  time  delay  on  the  Goff  drum  is  given  by  Goff  (see  Ref.  15, 
p.  20)  as  less  than  0.2  percent.  The  error  in  setting  the  zero  time  delay  for  the  ambiguity  dia¬ 
grams  is  about  50psec.  The  total  delay  is  24  msec,  so  the  error  introduced  by  setting  the  zero 
is  about  0.2  percent. 

4.4.5  Integrator  and  Sample-and-Hold  Error 

The  way  that  the  integrator  is  connected,  to  cancel  out  the  DC  drift  of  the  multiplier,  intro¬ 
duces  some  error  in  the  output.  A  small  amount  of  60-cps  noise  is  introduced  into  the  signal 
by  the  Goff  preamplifiers.  This  noise  is  squared  by  the  multiplier  and  results  in  a  small  amount 
of  120-cps  ripple  on  the  output  of  the  multiplier.  This  120-cps  ripple  offsets  the  zero  of  the 
integrator  slightly.  The  amount  of  error  in  the  output  due  to  this  offset  is  estimated  at  about 
1  percent. 

In  the  sample-and-hold  circuit,  with  the  sampling  relay  open,  a  slight  change  in  output  is 

encountered  as  the  grid  current  of  the  operational  amplifier  very  slowly  charges  or  discharges 
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the  sampling  capacitor.  Since  the  grid  current  is  on  the  order  of  10  amp  and  the  holding 
capacitor  is  a  0.1— jaf  capacitor,  the  output  voltage  changes  at  a  rate  of  about  1  mv/sec.  The  out¬ 
put  range  of  the  sample-and-hold  circuit  is  ±50  volts  and  the  hold  time  is  usually  l/4sec.  This 
drift  rate  corresponds  to  an  error  in  reading  of  0.00025  percent,  far  less  than  the  readability  of 
the  computer's  output  plotter. 

4.4.6  Summary  of  Errors 

The  errors  are  grouped  into  three  classes:  errors  producing  changes  in  the  height  of  the 
ambiguity  diagram,  errors  producing  changes  in  the  position  of  a  point  along  the  time  delay  axis, 
and  errors  producing  changes  in  the  position  of  a  particular  Doppler  shift  slice.  The  errors  are 
summarized  in  Table  III,  which  lists  the  maximum  errors  that  occur.  Since  the  contributions 
to  the  function  error  occur  on  both  channels,  the  individual  errors  are  counted  twice.  The  total 
estimated  error  in  the  height  of  the  ambiguity  diagram  is  5.4  percent  of  the  maximum  value; 
the  total  error  in  the  time  delay  is  0.54  percent  of  the  full-scale  value;  the  total  error  in  the 
position  of  a  particular  Doppler  shift  slice  is  4.4  percent  of  the  full-scale  value. 
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5.  AMBIGUITY  DIAGRAMS  OF  BAT  SIGNALS 


5.1  INTRODUCTION 

Signals  of  two  species  of  vespertilionid  bats,  Myotis  lucifugus  and  Lasiurus  borealis,  were 

obtained  while  they  were  in  pursuit  of  thrown  targets.  The  signals  of  the  Myotis  were  obtained 

in  the  laboratory  of  F.  A.  Webster;  the  signals  of  the  Lasiurus  were  obtained  in  the  field.  The 

Myotis  was  an  animal  that  had  been  in  captivity  for  several  months  and  had  been  trained  to  catch 

mealworms  thrown  into  the  air.  He  was  very  expert  at  catching  the  mealworms  and  missed  less 

than  10  percent  of  the  mealworms  thrown.  The  Lasiurus  was  an  animal  that  had  been  trained  by 

Dr.  A.  E.  Treat  and  his  family  to  come  close  to  their  house  to  receive  moths  that  were  thrown 

into  the  bat’s  approximate  path.  Since  some  of  the  moths  are  able  to  detect  the  ultrasonic  cries 
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of  the  bats  and  subsequently  to  avoid  the  bat,  the  Lasiurus  had  a  much  more  difficult  pursuit 
problem  than  the  Myotis.  The  Lasiurus,  however,  was  still  able  to  catch  more  than  50  percent 
of  the  moths  that  were  thrown  to  him. 

In  connection  with  certain  research  projects,  attempts  were  being  made  to  obtain  photographs 
of  the  bats  in  the  act  of  pursuing  and  catching  targets.  In  early  attempts,  a  room  was  flooded 
with  fruit  flies  and  bats  were  allowed  to  fly  about  and  catch  as  many  flies  as  they  could.  Since 
a  bat  is  small  and  travels  at  a  comparatively  high  speed,  most  of  the  pictures  that  were  obtained 
were  either  blank  (the  bat  was  out  of  the  field)  or  were  out  of  focus.  Even  when  a  sharp  photo¬ 
graph  of  the  bat  was  obtained,  the  bat  was  most  likely  just  flying  around  and  was  nowhere  near 
a  fruit  fly.  Because  of  the  limitations  imposed  by  the  size  of  the  photographic  field  and  the  depth 
of  field  of  photographic  equipment,  drastic  measures  had  to  be  taken  to  obtain  any  pictures  of 
bats  in  the  act  of  catching  targets. 

Since  the  bat  ate  mealworms  (the  larvae  of  the  mealworm  beetle)  out  of  a  dish  to  supplement 
what  was  caught  on  the  wing,  it  was  decided  that  an  attempt  would  be  made  to  teach  the  bat  to 
catch  thrown  mealworms,  which  could  be  more  easily  controlled  than  the  fruit  flies.  About  one 
bat  in  2  5  finally  learned  to  catch  mealworms  that  were  tossed  to  him  while  he  was  flying.  After 
learning  to  catch  these  targets,  the  bat  would  assume  a  regular  flight  pattern  around  the  area  of 
the  room  where  the  mealworms  were  thrown.  Cameras  could  be  set  up  to  photograph  the  bat 
since  his  exact  position  could  be  predicted  from  the  trajectory  of  the  thrown  mealworm.  Finally, 
a  ’’mealworm  mortar”  was  designed  and  built  which  fired  a  mealworm  to  an  exact  location  at  an 

exact  time.  This  mortar  made  possible  high-quality  photographs  of  the  bats  catching  targets. 
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Subsequently,  high-speed  motion  pictures  with  synchronized  sound  tracks  were  obtained. 

Thrown  moths,  however,  are  not  as  cooperative  as  thrown  mealworms.  When  the  moths 
are  thrown  into  the  air,  they  try  to  evade  the  bat  rather  than  going  to  the  position  in  the  photo¬ 
graphic  field  for  which  they  were  intended.  This  makes  photography  of  bats  catching  moths  very 
difficult.  The  bat,  however,  is  still  forced  to  fly  into  a  particular  area  in  order  to  catch  the  • 
thrown  moth;  consequently,  it  is  a  simple  matter  to  obtain  tape  recordings  of  his  ultrasonic  cries. 

The  bat's  echolocating  cries  were  picked  up  with  a  special  ultrasonic  microphone  that  was 

built  by  Dr.  L.  P.  Granath  of  Worcester  Polytechnic  Institute.  The  microphone  is  similar  in 

design  to  the  solid-dielectric  condenser  microphones  described  by  Kuhl,  Schodder,  and 
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Schroder.^  The  particular  microphones  used  have  a  frequency  response  that  extends  up  to 
200  kcps.  The  signals  from  the  microphone  are  recorded  on  a  high-speed  instrumentation  tape 
recorder  (Ampex  FR-1100)  running  at  60ips.  The  tape  recorder  has  a  frequency  response  that 
is  fairly  flat  out  to  lOOkcps  and  then  drops  very  rapidly. 


Since  the  microphone  that  is  used  is  not  perfect,  it  will  change  the  nature  of  the  bat  signal 
slightly.  Off  the  main  axis  of  the  microphone,  the  microphone  has  zeros  in  its  frequency  re¬ 
sponse.  The  general  response  of  the  microphone  is  roughly  given  by 


S(w) 


2  J^(aci)  sin  G) 
acn  sin  0 


(5-1) 


where  J^(x)  is  the  first-order  Bessel  function.  Since  the  position  of  the  bat  during  pursuit  of  the 
target  is  known,  the  microphone  is  aimed  in  such  a  way  that  the  bat  flies  on  the  axis  of  the  micro¬ 
phone  most  of  the  time.  This  is  not  too  difficult,  since  the  beamwidth  of  the  microphone  at  the 
highest  frequency  that  the  bat  uses  is  on  the  order  of  30°. 


5.2  MYQTIS  LUCIFUGUS 
5.2.1  Introduction 

The  bat  Myotis  lucifugus  is  the  little  brown  bat  common  to  the  New  England  area.  The  full- 
grown  adult  has  a  maximum  wing  spread  of  about  10  inches,  and  weighs  between  5  and  12  grams. 
The  bat’s  diet  consists  largely  of  insects  that  are  caught  in  flight.  He  has  some  eyesight  but  this 
is  barely  more  than  a  discrimination  between  light  and  no  light.  If  he  is  placed  in  a  dish  contain¬ 
ing  food,  too  close  for  him  to  make  use  of  his  sonar,  he  moves  around  in  a  random  fashion  until 
his  nose  actually  bumps  into  some  of  the  food  before  he  realizes  where  it  is  in  the  dish. 

The  particular  Myotis  whose  signals  appear  in  this  thesis  is  a  bat  that  was  caught  in  the  wild 
and  then  trained  to  catch  thrown  targets  in  the  laboratory.  The  laboratory  is  a  quonset  hut  that 
was  converted  by  F.  A.  Webster  into  a  place  to  study  bats.  The  bat  was  trained  to  catch  thrown 
targets  by  repeatedly  throwing  mealworms  into  his  flight  path.  When  he  had  learned  that  the  tar¬ 
gets  were  edible,  he  learned  to  fly  in  a  particular  flight  pattern  that  brought  him  over  the  area 
of  the  room  where  the  targets  were  thrown.  He  would  continue  to  fly  over  this  area  until  his  food 
was  thrown  to  him. 

One  or  two  microphones,  placed  in  the  path  of  his  flight,  were  used  to  obtain  tape  recordings 
of  his  sonar  signals  while  he  was  in  pursuit  of  a  target.  Figure  31  illustrates  the  wide  variation 
in  the  type  of  signal  used  in  pursuit  of  his  target.  It  is  an  oscillographic  record  made  from  tape 
recordings.  The  upper  trace  is  the  sound  pressure  on  one  microphone,  the  lower  trace  is  the 
sound  pressure  on  a  second  microphone,  simultaneously  recorded.  The  distortions  of  the  ampli¬ 
tude  waveform  are  partially  due  to  the  gain  pattern  of  the  bat's  emitted  signal,  the  gain  pattern 
of  the  microphones,  and  inclusion  of  additive  multipath  echoes.  The  position  of  the  bat  can  be 
determined  relative  to  the  positions  of  the  two  microphones  by  measuring  the  relative  time  delay 
between  a  particular  pulse  on  the  two  signal  tracks. 

Five  representative  pulses  were  chosen  from  this  sequence,  on  the  basis  of  being  ’’typical” 
and  having  a  high  signal -to-noise  ratio.  These  pulses,  assigned  the  numbers  2,  3,  4,  5,  and  6, 
as  indicated  in  Fig.  3  5,  were  analyzed  in  detail  by  means  of  the  ambiguity  diagram  computer. 


5.2.2  Ambiguity  Diagrams 

The  first  signal  to  be  analyzed  is  designated  as  pulse  2  in  Fig.  31.  A  more  detailed  photo¬ 
graph  of  the  pulse  is  shown  in  Fig.  32  along  with  a  plot  of  the  instantaneous  period  of  the  pulse. 
This  ’’instantaneous  period”  is  the  time  between  positive -going  zero  crossings  of  the  sound 
pressure  waveform.  It  can  be  thought  of  as  the  inverse  of  the  instantaneous  frequency  of  the 
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Fig.  31.  Sounds  emitted  by  Myotis  lucifugus  while  cruising  in  laboratory  room,  then  in  pursuit  of  target,  as  picked  up  by  two 
microphones  and  recorded  on  continuous  oscillographic  record.  Vertical  I ines  are  time  markers,  spaced  1  msec  apart.  Length 
of  each  strip  is  70  msec,  except  one  that  was  cut  short  to  avoid  splitting  pulse. 


signal.  The  snow  effect  that  is  seen  in  this  display  when  a  pulse  is  not  present  is  due  to  the  fact 
that  the  period  meter  is  being  triggered  by  noise  with  random  zero  crossings. 

Pulse  2  is  almost  sinusoidal  with  a  period  that  is  steadily  increasing  in  a  slightly  nonlinear 
fashion.  This  pulse  is  designated  as  a  cruising  pulse.  The  bat  uses  a  pulse  like  this  when  he  is 
cruising  in  search  of  food,  but  has  not  yet  located  any  targets.  The  starting  frequency  of  this 
pulse  is  lOOkcps,  the  final  frequency  is  40kcps,  and  the  duration  is  3.2  msec. 

A  three-dimensional  drawing  of  the  ambiguity  diagram  for  this  pulse  is  shown  in  Fig.  33. 

The  width  of  the  area  of  ambiguity,  if  the  velocity  of  the  target  is  known  to  be  zero,  is  about 
85|jLsec.  One  of  the  interesting  things  about  this  ambiguity  diagram  is  that  the  fine  structure  of 
the  ambiguity  surface  is  on  approximately  the  same  scale  as  the  coarse  structure.  If  the  signal- 
to-noise  ratio  were  high  enough,  and  the  velocity  of  the  target  were  known,  the  measurement 
accuracy  of  the  target  range  would  be  limited  only  by  the  width,  in  the  time -delay  direction,  of 
the  central  spike  of  this  surface.  If  the  same  modulation  were  used  on  a  much  higher  frequency 
carrier,  the  area  of  ambiguity  would  be  determined  by  roughly  the  upper  envelope  of  this  surface 
of  ambiguity.  The  upper  envelope  of  this  surface  can  be  visualized  as  the  surface  that  would 
result  if  a  silk  scarf  were  laid  over  the  entire  surface.  The  ratio  of  the  width  of  the  envelope 
to  the  width  of  the  central  spike  for  this  particular  ambiguity  diagram  is  approximately  7:1. 

The  quantity  l//3  is  equal  to  —5.3  msec.  This  means  that  if  the  acceleration  of  the  target  is 
small,  the  position  that  the  target  had  5.3  msec  before  the  echo  hit  it  is  known  as  accurately  as 
if  the  target  were  stationary.  This  is  not  quite  true,  because  the  signal  function  exhibits  a  curi¬ 
ous  spreading  behavior  when  the  target  velocity  is  on  the  order  of  10  m/sec.  The  3-db  width  of 
the  signal  function  envelope  for  a  target  velocity  of  lOm/sec  is  about  llOpsec.  This  width  cor¬ 
responds  to  a  range  accuracy  of  1.9  cm. 

The  second  signal  that  is  analyzed  is  shown  in  Fig.  34  and  is  designated  as  pulse  3.  The 
position  of  this  signal  in  the  over-all  pulse  train  can  be  seen  in  Fig.  31.  This  pulse  occurs 
either  just  before  or  just  after  the  bat  detects  the  target.  Detection  of  the  target  is  followed  by 
the  decrease  in  the  interpulse  period  that  can  be  seen  in  Fig.  31  and  is  plotted  in  Fig.  3  5.  The 
interpulse  period  varies  during  cruising  between  20  and  100  msec.  When  the  target  is  detected, 
the  interpulse  period  drops,  the  duration  of  the  bat's  signal  decreases,  and  the  rate  of  change 
of  frequency  increases.  In  this  particular  plot,  the  interpulse  period  is  relatively  constant  and 
then  decreases  sharply.  The  constant  interpulse  period  before  detection  of  the  target  is  not 
typical  of  the  general  case.  Usually,  the  interpulse  period  varies  up  and  down  in  a  rather  irregu¬ 
lar  fashion  until  the  target  is  detected;  the  interpulse  period  then  decreases  as  is  shown  in  Fig.  35. 
The  jump  in  the  interpulse  period  just  before  the  long  sequence  of  closely  spaced  pulses  is  pres¬ 
ent  in  many  cases.  It  has  been  hypothesized  that  this  is  due  to  the  bat's  taking  in  a  breath  of  air 
so  that  he  can  emit  the  entire  sequence  of  closely  spaced  pulses  with  one  breath.  In  this  case, 
the  relatively  constant  amount  of  interpulse  period  before  detection  might  be  due  to  the  fact  that 
the  bat  knew  that  a  mealworm  was  being  thrown  and  was  searching  a  known  area  for  the  target. 

A  three-dimensional  drawing  of  the  ambiguity  diagram  for  pulse  3  is  shown  in  Fig.  36.  The 
over-all  nature  of  this  surface  is  very  similar  to  the  ambiguity  diagram  for  pulse  2.  This  is 
interesting  since  the  shape,  the  duration,  and  the  frequency  modulation  of  pulse  3  are  somewhat 
different  from  those  of  pulse  2.  A  close  examination  of  the  structure  of  the  ambiguity  diagram 
for  pulse  3  reveals  that  the  parameters  of  this  surface  are  different.  The  width  of  the  v  =  0  slice 
of  the  surface  is  70psec;  the  width  of  the  v  =  lOm/sec  slice  is  100  psec;  the  quantity  l//3  is  equal 
to  —3.5  msec. 
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Fig.  33.  Cruising-pulse  ambiguity  diagram  for  Myotis  lucifugus  (pulse  2). 


Fig.  34.  Myotis  lucifugus  locating 
pulse  (pulse  3). 
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Fig.  35.  Time  interval  between  a  pulse  and  preceding  pulse  for  record 
shown  in  Fig.  31.  Small  numbers  inside  circles  correspond  to  numbers 
assigned  to  individual  pulses. 
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Fig.  36.  Locating-pulse  ambiguity  diagram  for  Myotis  lucifugys  (pulse  3). 
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Fig.  37.  Myotis  lucifugus  pursuit  pulse  (pulse  4). 

Pulse  4  is  emitted  after  the  bat  has  located  the  target.  A  photograph  of  pulse  4  is  shown 
in  Fig.  37.  The  upper  frequency  of  pulse  4  is  about  one -half  of  what  it  was  for  pulses  2  and  3, 
and  the  duration  has  decreased  from  3.2  msec  in  pulse  2  and  2.2  msec  in  pulse  3  to  0.8  msec  in 
pulse  4.  A  three-dimensional  drawing  of  the  ambiguity  diagram  for  pulse  4  is  illustrated  in 
Fig.  38.  The  width  of  the  v  =  0  slice  of  the  surface  is  85  psec;  the  width  of  the  v  =  10  m/sec 
slice  is  110  psec;  the  quantity  l//3  is  equal  to  —1.6  msec.  This  ambiguity  diagram  illustrates 
the  general  trend  that,  as  the  bat  approaches  the  target,  he  uses  a  pulse  that  gives  him  a  smaller 
amount  of  range  error  for  a  given  velocity  error.  Although  the  ambiguity  diagram  for  pulse  2 
has  nearly  30  maxima  (for  the  amount  of  diagram  illustrated)  and  the  ambiguity  diagram  for 
pulse  3  has  15  maxima,  the  ambiguity  diagram  for  pulse  4  has  only  7  maxima. 

Pulses  5  and  6  are  emitted  when  the  bat  is  very  close  to  the  target.  He  knows  where  the 
target  is  and  is  keeping  tabs  on  it  to  make  sure  that  it  doesn't  change  course.  These  two  pulses 
are  illustrated  in  Fig.  39(a-b),  respectively.  The  instantaneous  frequency  of  these  pulses  is  no 
longer  defined  in  the  relatively  simple  way  that  the  instantaneous  frequency  of  the  three  former 
pulses  has  been  defined.  The  signals  are  no  longer  sinusoidal  in  shape,  but  contain  a  large 
amount  of  second  harmonic.  If  the  instantaneous  frequency  of  these  signals  is  defined  as  the 
inverse  of  the  time  between  the  zero  crossings  of  the  fundamental  waveform,  upper  and  lower 
frequency  bounds  can  be  defined.  The  frequency  range  for  pulse  5  is  40  kcps  down  to  20 kcps; 
the  duration  is  0.4  msec.  The  frequency  range  for  the  fundamental  of  pulse  6  is  50  kcps  down  to 
17  kcps;  the  duration  is  0.3  msec.  Three-dimensional  drawings  of  the  ambiguity  diagrams  for 
pulses  5  and  6  are  illustrated  in  Fig.  40(a-b),  respectively. 

The  ambiguity  diagram  of  pulse  5  has  only  3  maxima.  The  width  of  the  v  =  0  slice  of  the 
surface  is  lOOpsec;  the  width  of  the  envelope  of  the  v  =  10  m/sec  slice  of  the  surface  is  100  psec; 
the  quantity  1  //3  is  equal  to  —0.5  msec.  The  ambiguity  diagram  of  pulse  6  has  only  one  maximum. 
The  width  of  the  v  =  0  slice  of  this  surface  is  125  psec;  the  width  of  the  v  10  m/sec  slice  of  this 
surface  is  12  5  psec;  the  quantity  \/p  is  equal  to  —0.4  msec. 

5.2.3  Discussion 

A  summary  of  the  parameters  of  the  pulses  and  of  the  ambiguity  diagrams  is  contained  in 
Table  IV.  It  is  interesting  that  the  one -Napier  width  of  the  envelope  of  the  waveform  in  the 
v  =  0  slice  of  the  ambiguity  diagram  remains  relatively  constant  for  the  various  pulses.  This 
means  that  the  amount  of  resolution  that  the  bat  has  is  relatively  independent  of  the  particular 
signal  that  he  uses,  out  of  the  class  of  signals  that  are  available  to  him. 
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Fig.  38.  Pursuit-pulse  ambiguity  diagram  for  Myotis  lucifugus  (pulse  4). 
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Fig.  39.  Myotis  lucifugus  pursuit  pulse. 
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Fig.  40.  Pursuit-pulse  ambiguity  diagram  for  Myotis  lucifugus. 
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TABLE  IV 

SUMMARY  OF  THE  PARAMETERS  OF  BAT  PULSES  AND  AMBIGUITY 

DIAGRAMS 

BAT  PULSE 

AMBIGUITY  DIAGRAM 

Duration 

Start  f 

Stop  f 

v  =  0  width 

v  =  10  width 

i/p 

Pulse  Number 

(msec) 

(kcps) 

(kcps) 

(psec) 

(psec) 

(msec) 

Maxima 

Myotis  lucifugus 

2 

3.2 

100 

40 

85 

110 

-5.3 

30 

3 

2.2 

100 

35 

70 

100 

-3.5 

15 

4 

0.8 

50 

30 

85 

110 

-1.6 

7 

5 

0.4 

40 

22 

100 

100 

-0.5 

3 

6 

0.3 

(50) 

(17) 

125 

125 

-0.4 

1 

Lasiurus  borealis 

7 

9.6 

60 

35 

380 

500 

-70 

300 

Another  interesting  observation  can  be  made  about  how  the  quantity  i/p  behaves.  For  the 
cruising  pulse,  pulse  2,  i/p  is  equal  to  —5.3  msec.  After  the  bat  detects  the  target  and  closes 
in  on  it,  the  magnitude  of  i/p  decreases  until  i/p  is  less  than  —0.5msec  for  pulses  5  and  6. 

The  meaning  that  was  associated  with  the  quantity  i/p  in  Sec.  2.3.3  is  that  this  is  the  time  for 
which  the  position  of  the  target  is  determined,  relative  to  the  time  that  the  echo  bounces  off  the 
target. 

If  it  is  hypothesized  that  the  magnitude  of  i/p  is  equal  to  the  time  that  it  takes  for  the  bat’s 
emitted  pulse  to  travel  through  the  air  to  the  target,  it  then  follows  that  the  time  at  which  the 
pulse  is  emitted  is  the  time  for  which  the  position  of  the  target  is  determined.  This  means  that 
the  position  of  the  target,  for  the  time  at  which  the  pulse  is  emitted,  is  given  at  the  output  of 
the  ideal  receiver,  independent  of  the  velocity  of  the  target. 

In  order  to  check  this  hypothesis,  the  values  given  in  Table  IV  for  i/p  can  be  examined  to 
determine  if  they  correspond  to  a  reasonable  bat-to-target  distance.  Multiplying  the  values  of 
i/p  by  the  velocity  of  sound  in  air,  346  m/sec,  the  distance  for  the  cruising  pulse  is  1.8  in;  for 
the  detecting  pulse,  1.2  m;  and  for  the  last  pursuit  pulse,  it  is  0.14  m.  These  values  seem  to 
be  a  little  high,  but  further  work  is  needed  to  determine  whether  or  not  they  are  actually  incon¬ 
sistent  with  the  above  hypothesis. 

Since  the  power  of  the  echo  decreases  as  the  inverse  fourth  power  of  the  distance  from  the 
bat  to  the  target,  it  seems  logical  that  the  bat  would  be  willing  to  accept  some  loss  in  range 
resolution  during  cruising  in  order  to  be  able  to  use  a  more  powerful  signal.  Suppose,  for 
example,  that  the  bat  can  detect  a  target  at  a  distance  of  1  m.  If  a  pulse  that  has  the  same 

A 

energy  is  used,  when  he  is  0.1  m  from  the  target,  the  echo  energy  is  10  greater! 

For  the  particular  pulse  train  that  has  been  analyzed,  the  duration  of  the  cruising  pulse  is 
more  than  ten  times  the  duration  of  the  final  pursuit  pulse  and  the  amplitude  of  the  cruising 
pulse  is  about  three  times  the  amplitude  of  the  pursuit  pulse.  Using  these  figures,  the  energy 
contained  in  the  cruising  pulse  is  100  times  the  energy  contained  in  the  pursuit  pulse.  This 
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means  that  if  the  bat  had  a  signal-energy-to-noise-power-density  ratio  of  one  during  cruising, 

2 

he  would  have  a  ratio  of  10  when  he  is  on  the  final  leg  of  the  pursuit.  Referring  back  to  the 
ambiguity  diagram  for  pulse  6,  Fig.  40(b),  the  type  of  pulse  used  in  the  final  leg  of  the  pursuit, 
it  can  be  seen  that  with  this  very  high  signal-to-noise  ratio,  the  position  of  the  target  is  deter¬ 
mined  only  by  the  very  top  of  the  central  spike  of  this  surface.  This  should  lead  to  extremely 
good  range  accuracy. 

5.3  LASIURUS  BOREALIS 

5.3.1  Introduction 

The  bat  Lasiurus  borealis,  commonly  known  as  the  red  bat,  grows  to  have  a  wingspread  of 
about  18  inches,  and  weighs  between  10  and  30  grams. 

The  particular  bat  whose  signals  appear  in  this  report  was  trained  by  Dr.  A.  E.  Treat  and 
his  family  to  catch  moths  thrown  in  the  air.  The  tape  recordings  were  made  at  Dr.  Treat’s 
summer  home,  at  Tyringham,  Massachusetts,  during  an  evening  when  the  bat  was  being  thrown 
moths  by  several  different  people. 

A  microphone  was  situated  in  a  fixed  position  near  the  house,  pointing  in  the  general  direc¬ 
tion  of  the  bat’s  flight  path.  The  same  troubles  with  signal-to-noise  ratio  that  were  encountered 
with  the  signals  from  the  Myotis  were  encountered  with  the  signals  from  the  Lasiurus,  with  the 
exception  that  since  the  signals  were  recorded  in  the  free-space  outdoors,  there  were  no 
multipath-echo  problems.  The  signal  shown  was  chosen  on  the  basis  of  having  a  good  signal-to- 
noise  ratio  and  is  typical  of  the  Lasiurus  pulses. 

5.3.2  Ambiguity  Diagrams 

The  signal  that  is  analyzed  is  shown  in  Fig.  41  and  is  designated  as  pulse  7.  This  is  a 
cruising  pulse.  That  is,  it  was  emitted  while  the  Lasiurus  was  in  search  of  a  target  but  had 
not  yet  detected  it.  This  Lasiurus  pulse  is  very  different  from  any  of  the  Myotis  pulses.  The 
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Fig.  41.  Lasiurus  borealis  cruising  pulse  (pulse  7). 

frequency  starts  at  about  60  and  sweeps  to  35kcps.  The  sweep  is  very  nonlinear  and  exhibits 
an  almost-constant-frequency  portion  at  the  end  of  the  pulse.  The  pulse  has  an  extremely  long 
duration,  9.6  msec,  more  than  twice  the  duration  of  the  longest  Myotis  pulse.  The  autocorrela¬ 
tion  function  for  this  signal  is  shown  in  Fig.  42.  Only  half  the  function  is  plotted,  since  the 
function  is  symmetrical  around  t  =  0.  This  autocorrelation  function  has  more  than  four  times 
the  width  of  the  autocorrelation  function  of  the  Myotis  cruising  pulse.  The  large  number  of 
peaks  and  valleys  seem  to  indicate  that  this  bat  cannot  make  use  of  the  fine  structure  information 
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indicated  in  this  autocorrelation  function.  This  is  not  an  exact  statement,  since,  with  a  high 
enough  signal-to-noise  ratio,  the  position  of  the  target  will  be  determined  by  the  position  of  the 
very  highest ‘peak  of  the  ambiguity"  diagram.  A  drawing  of  the  ambiguity  diagram  is  illustrated 
in  Fig. 43 . 

Since  the  fine  structure  differs  so  greatly  from  that  of  the  Myotis  slices,  a  different  way  of 
representing  the  ambiguity  diagram  for  pulse  7  was  devised.  For  the  Myotis,  the  entire  ambi¬ 
guity  surface  was  shown,  represented  as  an  opaque  surface,  with  the  connections  between  the 
individual  slices  filled  in. 

For  the  Lasiurus,  however,  only  the  envelopes  of  the  individual  slices  are  represented, 
with  the  area  inside  these  envelopes  filled  in  with  vertical  lines  representing  the  oscillations  of 
the  waveforms  of  the  individual  slices.  The  connections  between  the  slices  have  been  omitted 
because  of  the  difficulty  in  drawing  these  lines  and  because  there  really  isn't  enough  information 
to  do  an  intelligent  job  of  filling  in  the  blank  volume.  Drawn  in  this  way,  some  of  the  lower  sur¬ 
face  is  also  revealed,  and  one  may  observe  the  shape  of  the  small  sections  that  split  off  beyond 
the  Doppler-equal-to-zero  slice,  which  would  otherwise  be  hidden  by  the  "mountain  range"  of 
peaks. 

The  distribution  of  the  area  of  ambiguity  is  so  different  from  that  of  the  Myotis,  that  differ¬ 
ent  scales  are  used  on  both  the  time -delay  and  the  Doppler  axis  of  the  ambiguity  diagram.  The 
relatively  restricted  range  for  the  Doppler  scale,  —5  to  +5  m/sec,  was  used  because  of  the  time- 
delay  limitations  of  the  ambiguity  diagram  computer  and  not  because  the  Lasiurus  flies  slower 
or  catches  slower  targets  than  the  Myotis  does.  On  the  contrary,  the  Lasiurus  usually  flies 
about  twice  as  fast  as  the  Myotis. 

Rather  than  exhibiting  the  relatively  smooth  behavior  that  the  Myotis  ambiguity  diagrams 
exhibited,  this  ambiguity  diagram  splits  into  several  sections  off  the  Doppler-equal-to-zero 
slice.  This  is  probably  due  to  the  fact  that  the  signal  does  not  have  linear  frequency  modulation 
but  has  a  very  nonlinear  type  of  frequency  modulation.  The  ambiguity  diagram  consists  of  the 
sum  of  the  contributions  of  the  various  parts  of  the  signal,  as  indicated  by  Eq.  (2-27).  When  the 
signal  is  shifted  from  the  r  =  0  point,  each  portion  of  the  signal  is  multiplied  by  a  portion  of  the 
signal  in  its  near  vicinity.  Therefore,  the  portion  of  the  signal  in  the  near  vicinity  of  any  por¬ 
tion  is  the  only  part  that  has  a  major  effect  on  this  portion's  contribution  to  the  area  of  ambiguity. 

The  width  of  the  zero  Doppler  slice  of  this  surface  is  380|jLsec;  the  width  of  the  — 5-m/sec 
slice  is  500  psec  (the  l/e  width);  the  quantity  l//3  is  equal  to  70  msec.  The  range  accuracy  cor¬ 
responding  to  a  time-delay  accuracy  of  380  psec  is  6.6  cm.  Since  this  long  pulse  is  used  only 
for  long-range  detection,  it  is  not  inconsistent  that  the  bat  does  not  have  the  excellent  range 
measurement  accuracy  that  the  Myotis  had. 

The  nature  of  the  Lasiurus  pulses  changes  greatly  as  the  bat  closes  in  on  his  target.  Simi¬ 
lar  to  the  Myotis,  the  Lasiurus  shortens  his  pulses  until  they  are  only  1  msec  in  duration  when 
he  is  in  the  final  leg  of  his  pursuit.  The  Lasiurus,  however,  does  not  decrease  the  frequency 
of  his  signals  as  the  Myotis  does.  He  sweeps  over  approximately  the  same  frequency  range  in 
the  shorter  length  of  time.  The  pursuit  signals  of  the  Lasiurus  are  not  included  here  because 
ambiguity  diagrams  have  not  yet  been  computed  for  them. 


61 


Fig.  42. 


Cruising-pulse  autocorrelation  function  for  Lasiurus  boreal  is  (pulse  7). 
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Fig.  43.  Slices  of  cruising-pulse  ambiguity  diagram  for  Lasiurus  boreal  is  (pulse  7). 
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6.  CONCLUSIONS 


6.1  SUMMARY 


It  has  been  the  purpose  of  this  study  to  determine  the  range  and  range-rate  measurement 
capabilities  of  various  bat  sonar  signals.  The  results  indicate  that: 

(a)  It  is  possible  to  derive  a  wide-band  echolocation  theory  that  predicts, 

on  the  basis  of  the  a  posteriori  model,  the  range  and  range-rate  accuracy 
of  a  given  signal. 

(b)  Accurate  angular  measurements  can  be  made  with  the  use  of  wide-band 
signals  and  gain-dispersive  antennas. 

(c)  An  analog  ambiguity  diagram  computer  that  calculates  the  signal 
function  characterizing  the  operation  of  the  a  posteriori  model  can 
be  built  and  operated  satisfactorily. 

(d)  The  effective  area  of  ambiguity  of  this  signal  function  is  dependent  on 
the  bandwidth  of  the  transmitted  signal  and  is  given  approximately  by  the 
relation 
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(e)  The  signals  that  are  used  by  the  bats  Mvotis  lucifugus  and  Lasiurus 
borealis .  pulse  by  pulse,  have  almost  no  Doppler  information,  but 
have  extremely  good  range-resolution  capabilities. 


Although  bats  may  not  use  the  principles  of  the  a  posteriori  model  in  a  manner  directly 
analogous  to  the  one  presented,  it  can  be  stated  that  a  bat  cannot  do  better  than  the  "ideal 
receiver"  under  the  conditions  of  additive  Gaussian  noise  that  are  stated  in  Chapter  2.  This 
derivation  places  a  limit  on  the  capabilities  of  any  man-made  or  biological  echolocation  system. 

In  particular,  this  ideal-receiver  approach  demonstrates  that  the  signals  used  by  Myotis 
and  Lasiurus  are  well  suited  for  pulse  compression.  The  3.2-msec-long  cruising  pulse  of  the 
Myotis,  corresponding  to  a  pulse  length  in  air  of  1.1m,  compresses  when  passed  through  the 
ideal  receiver  to  a  width  of  85psec,  which  in  turn  corresponds  to  a  resolution  of  1.5  cm. 

The  nature  of  the  ambiguity  diagrams  that  were  computed  for  the  Myotis  pulses  demonstrate 
that  this  bat  cannot  do  much,  if  any,  Doppler  discrimination  on  a  single-pulse  basis.  The 
height  of  the  surface  of  ambiguity  is  nearly  as  high  for  the  maximum  a  priori  velocity  as  it  is 
for  zero  velocity.  This  indicates  that  velocity  determination  on  a  single-pulse  basis  is  only 
possible  under  very  high  signal-to-noise-ratio  conditions.  However,  the  bat  does  not  use  only 
a  single  pulse  but  in  fact  uses  a  train  of  pulses.  If  on  each  pulse  he  has  good  range  information, 
the  train  of  pulses  will  give  him  both  range  and  velocity  information. 

On  the  other  hand,  the  ambiguity  diagram  computed  for  the  Lasiurus  cruising  pulse  demon¬ 
strates  that  this  bat  has  some  range  and  some  velocity  resolution.  The  380-fisec  width  of  the 
ambiguity  diagram  corresponds  to  a  range  resolution  of  6.5  cm.  The  peak  of  the  target-velocity- 
equal-to-zero  slice  is  substantially  higher  than  any  of  the  other  peaks  of  this  ambiguity  diagram, 
indicating  that  under  favorable  signal-to-noise-ratio  conditions,  the  bat  signal  is  capable  of  being 

used  to  measure  velocity.  The  ambiguity  diagram  for  this  particular  signal  is  similar  to  the 

22 

circularly  symmetrical  ambiguity  diagram  proposed  by  Klauder,  in  that  this  ambiguity  dia¬ 
gram  breaks  into  several  minor  ridges  off  the  target-velocity-equal-to-zero  slice. 

A  complete  treatment  of  the  bat  sonar  problem  would  include  the  treatment  of  the  problem 
of  simultaneously  determining  range,  range  rate,  azimuth,  elevation,  angular  velocity,  and 
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type  of  target.  Since  this  problem  seems  to  be  very  complex,  if  not  impossible  to  handle,  the 
simpler  problem  of  simultaneously  determining  range  and  elevation  is  treated.  The  results  of 
this  treatment  indicate  that  it  is  possible  to  determine  the  elevation  of  a  target  from  the  spec¬ 
trum  of  the  received  signal  if  a  gain-dispersive  antenna  is  used  as  the  receiving  aperture.  The 
two  simple  examples  treated  serve  to  illustrate  the  technique  and  the  possibilities  of  this  type  of 
system.  The  results  show  that  signals  with  equivalent  power  density  spectrum  give  identical 
performance.  As  in  all  a  posteriori  models,  the  signal-to-noise  ratio  affects  the  operation  of 
the  system  greatly. 

An  analog  ambiguity  diagram  computer  was  designed  and  built  to  calculate  the  signal  function 
defined  in  Chapter  2.  The  maximum  error  introduced  by  this  computer  is:  height,  5.4  percent; 
time  delay,  0.54  percent;  position  of  a  Doppler  shift  slice,  4.4  percent. 

6.2  SUGGESTIONS  FOR  FURTHER  INVESTIGATION 
6.2.1  Equipment  Improvements 

The  weakest  link  in  the  ambiguity  diagram  computer  is  the  input  tape  reproducer.  The 
amount  of  flutter  introduced  by  this  machine  limits  the  resolution  in  the  Doppler  direction  and 
restricts  the  use  of  the  ambiguity  diagram  computer.  It  is  doubtful  if  the  addition  of  a  servo 
system  on  the  present  tape  reproducer  could  reduce  the  flutter  by  an  appreciable  amount.  Since 
the  operation  of  this  particular  tape  reproducer  is  already  commensurate  with  the  state  of  the 
art  of  commercial  tape  reproducers,  it  is  not  reasonable  to  assume  that  a  different  commerical 
tape  reproducer  would  operate  significantly  better.  One  possible  solution  would  be  to  record 
the  signal  onto  another  magnetic  drum  and  then  re-record  it  onto  the  Goff  drum  with  the  Doppler 
shift.  Another  solution  would  be  to  record  the  Doppler- shifted  signal  onto  the  Goff  drum  by 
recording  the  signal  with  the  movable  head,  moving  at  a  relatively  high  rate,  while  at  the  same 
time  the  non- Doppler  -  shifted  signal  is  being  recorded  with  a  stationary  head.  A  special  drive 
mechanism  would  have  to  be  designed  in  order  to  move  the  movable  head  at  the  high  rates  nec¬ 
essary  to  provide  the  Doppler  shift.  By  recording  the  two  signals  at  the  same  time,  errors 
introduced  by  flutter  would  cancel  each  other  out. 

A  digital  computer  could  be  used  to  calculate  the  ambiguity  diagrams.  This  would  introduce 
a  versatility  not  enjoyed  by  the  analog  computer.  Using  the  presently  available  digital  computers, 
the  accuracy  of  the  ambiguity  diagram  would  probably  be  much  better  than  the  accuracy  obtained 
with  the  analog  computer.  The  entire  error  would  probably  be  introduced  by  the  first  conversion 
from  the  analog  to  the  digital  form  of  input  signal.  Calculation  of  an  ambiguity  diagram  by  a 
digital  computer  is  a  very  time-consuming  operation,  however,  and  would  probably  be  much 
more  expensive  than  with  the  analog  computer. 

The  60-cps  hum  in  the  Goff  magnetic-drum  reproduce  amplifiers  affects  the  operation  of 
zero  setting  the  output  of  the  multiplier  as  described  in  Chapter  4.  Elimination  of  this  hum, 
while  the  multiplier  zero  is  being  set,  would  eliminate  this  problem.  The  120-cps  hum,  resulting 
from  the  squaring  of  the  60-cps  hum  signals,  can  be  eliminated  by  grounding  the  two  multiplier 
inputs  during  the  time  that  the  zero  is  being  set.  These  inputs  could  very  easily  be  grounded 
with  relays  connected  to  the  integrate  one-shot  multiplier.  The  multiplier  inputs  would  then  be 
grounded  until  the  actual  integration  was  to  occur.  The  inputs,  of  course,  would  be  connected 
to  the  reproduce  amplifier  outputs  while  the  signal  is  being  multiplied  and  integrated. 

Some  signals,  for  which  it  is  desirable  to  obtain  ambiguity  diagrams,  cannot  be  analyzed 
with  this  ambiguity  diagram  computer  because  their  entire  ambiguity  diagram  lies  very  near 
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the  target-velocity-equal-to-zero  axis.  An  example  of  such  a  pulse  is  the  echolocation  pulse 
of  the  bat  Rhinolophus  ferrum-equinum.  This  pulse  has  a  long  constant-frequency  section,  as 
well  as  a  short  FM  section  at  the  ehd  of  the  pulse.  A  typical  center  frequency  for  this  pulse  is 
85kcps;  a  typical  duration  is  20msec.  With  this  signal,  the  entire  ambiguity  diagram  is  con¬ 
tained  between  plus  and  minus  0.1  m/sec,  way  out  of  the  range  of  this  ambiguity  diagram  com¬ 
puter.  Such  a  signal  could  be  analyzed  if  this  relatively  narrow-band  signal  were  played  at 
normal  speed  on  a  tape  reproducer,  beat  down  in  frequency  so  that  it  would  become  a  wide-band 
signal  in  another  frequency  range,  and  re-recorded  on  a  tape  recorder.  This  beating  down 
could  probably  be  accomplished  with  the  aid  of  an  oscillator,  the  multiplier  that  is  used  in  this 
computer,  and  a  band -pass  filter. 

6.2.2  Field  Measurements  of  the  Capabilities  of  Bat  Sonar  Systems 

This  research  project  has  built  up  the  theoretical  basis  for  the  comparison  of  the  perform¬ 
ance  of  the  sonar  system  of  individual  bats  to  the  performance  of  an  ideal  receiver.  Further 
measurements  of  the  capabilities  of  bats  are  needed  in  order  to  determine  if  their  performance 
is  indeed  comparable  to  that  of  the  ideal  receiver. 

Some  experiments  that  might  be  undertaken  are:  a  determination  of  the  resolution  of  targets 
in  the  presence  of  Gaussian  jamming  noise,  the  range  accuracy  of  the  bat's  system,  and  the 
angular  accuracy  of  the  bat's  system.  Such  experiments  might  be  relatively  difficult  to  devise 
and  carry  out. 

Simultaneous  high-speed  motion-picture  photography  and  tape  recording  can  be  used  to 
great  advantage  in  studying  the  bats.  As  an  example  of  its  application,  one  could  use  this 
technique  to  check  whether  or  not  the  quantity  c/j3  is  equal  to  the  distance  between  the  bat  and 
his  target  as  is  hypothesized  in  Chapter  5. 

6.2.3  Ambiguity  Diagrams  of  Various  Bat  Signals 

There  are  many  species  of  bats  that  use  many  different  echolocating  signals.  Each  bat  may 
have  a  special  reason  for  using  the  particular  signal  that  he  uses.  It  would  be  very  interesting 
to  know  what  kind  of  ambiguity  diagrams  result  from  the  pursuit  pulses  of  the  Lasiurus.  The 
fishing  bat,  Noctillio,  seems  to  switch  back  and  forth  between  CW  pulses  and  FM  pulses.  This 
bat  uses  these  signals  to  locate  fish  swimming  in  the  water.  The  long-eared  bats,  Plecotus 
townsendii,  emit  pulses  that  are  much  lower  in  frequency  than  the  signals  of  the  Myotis.  These 
signals  are  very  rich  in  harmonic  content  and  may  produce  a  very  different  type  of  ambiguity 
diagram.  The  number  of  different  species  of  bats  is  very  large.  It  would  be  interesting  indeed 
to  obtain  the  ambiguity  diagrams  of  all  of  them,  for  each  of  their  different  uses,  but  time  is  not 
available  to  pursue  all  of  them. 
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